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ABSTRACT 

We present a sensitive, multiwavelength submillimeter continuum survey of 153 
young stellar objects in the Taurus-Auriga star formation region. The submillime- 
ter detection rate is 61% to a completeness limit of ~10mJy (3-(t) at 850 //m. The 
inferred circumstellar disk masses are log-normally distributed with a mean mass of 
~ 5 X 10~^ M0 and a large dispersion (0.5 dex). Roughly one third of the submillimeter 
sources have disk masses larger than the minimal nebula from which the solar system 
formed. The median disk to star mass ratio is 0.5%. The empirical behavior of the sub- 
millimeter continuum is best described as oc iy2.o±o.5 bg^-^ggn 350 /xm and 1.3 mm, 
which we argue is due to the combined effects of the fraction of optically thick emission 
and a flatter frequency behavior of the opacity compared to the interstellar medium. 
This latter effect could be due to a substantial population of large dust grains, which 
presumably would have grown through collisional agglomeration. In this sample, the 
only stellar property that is correlated with the outer disk is the presence of a com- 
panion. We find evidence for significant decreases in submillimeter flux densities, disk 
masses, and submillimeter continuum slopes along the canonical infrared spectral en- 
ergy distribution evolution sequence for young stellar objects. The fraction of objects 
detected in the submillimeter is essentially identical to the fraction with excess near- 
infrared emission, suggesting that dust in the inner and outer disk are removed nearly 
simultaneously. 

Subject headings: circumstellar matter — planetary systems: protoplanetary disks — 
solar system: formation — stars: pre-main— sequence 



1. Introduction 

The formation and early evolution of stars are intimately coupled to the properties of their 
accompanying circumstellar disks of gas and dust. These disks also provide the material reservoirs 
for the assembly of planetary systems. Angular momentum conservation dictates that a collapsing 
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molecular cloud core with some initial rotation will result in both a central protostar and a flattened 
circumstellar disk (e.g., Terebey, Shu, &: Cassen 1984). Indirect observations indicate that disks are 
essentially ubiquitous in young star clusters, while optical images in silhouette (O'Dell & Wen 1994) 
and millimeter spectral line confirmations of Keplerian rotation (e.g., Simon, Dutrey, &: Guilloteau 
2001) provide more direct evidence in specific cases. Comparisons of infrared observations with 
physical models of young stellar objects (YSOs; here taken to mean a young star and its associated 
circumstellar material) have led to a sequence of evolutionary stages which occur before the start 
of the main-sequence (Lada & Wilking 1984; Adams & Shu 1986; Adams, Lada, & Shu 1987). In 
the Class I stage, an extended circumstellar envelope is rapidly dumping material onto a central 
protostar and a massive accretion disk. After the supply of envelope material is dissipated, the 
YSO becomes a Class II object, with a disk that is actively accreting material onto a central, 
optically visible star. In the final Class III stage, at least the inner part of the circumstellar 
disk has been evacuated, although the dominant physical mechanism for this process remains in 
debate (sec HoUcnbach, Yorke, & Johnstone 2000). The most interesting possibility, at least from a 
cosmogonical viewpoint, is that the gas and dust in the disk have agglomerated into larger objects 
in a developing planetary system. 

Observations of the morphology of the broadband spectral energy distribution (SED) and 
various diagnostics of accretion can be used to trace the evolution of a YSO. Longward of ~1 /xm, 
the SED of a YSO is composed of a continuum of thermal spectra from the radially distributed 
circumstellar dust, modified by the radiative transfer properties of the grains. Changes in the SED 
through the evolutionary sequence are indicative of the loss of circumstellar components in the 
system; first the envelope and then the disk. The slope of the infrared SED is determined by the 
radial temperature distribution of the circumstellar dust (e.g., Adams, Lada, &: Shu 1987; Beckwith 
et al. 1990). Therefore, measurements of infrared colors provide a relatively simple observational 
constraint on the temperature structure of a disk. However, more detailed physical interpretations 
of the infrared SED are challenging, due to the strong dependence on the relatively unknown 
radiative transfer properties of the grains and detailed disk structure (e.g., the inner disk radius or 
vertical scale height). 

In the early evolution stages (Class I and II), material from the inner disk is dragged in mag- 
netosphcric funnel flows to the stellar surface, with an accretion shock resulting upon impact (see 
the review by Najita et al. 2000). This process is responsible for the observed continuum excesses 
(Calvet & Gullbring 1998; Johns-KruU & Valenti 2001; Muzerolle et al. 2003) and the shapes and 
strengths of emission lines in YSOs (Hartmann, Hewett, & Calvet 1994; Muzerolle, Hartmann, & 
Calvet 1998; Muzerolle, Calvet, &; Hartmann 2001). The most common observational measurement 
providing a breakdown of objects as accreting or non- accreting is the equivalent width (W) of the 
Ha emission line. Although a standard division at = 10 A was set by historical instrument limi- 
tations rather than a physical motivation, this criterion provides an cfl"ectivc discriminant as many 
properties of weak-line (WTTSs; < lOA; non-accreting) and classical (CTTSs; W > 10 A; ac- 
creting) T Tauri stars are remarkably different (e.g., Ghez, Neugebauer, &; Matthews 1993; Osterloh 
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& Beckwith 1995; Chiang, Phillips, & Lonsdale 1996; Stelzer & Neuhauser 2001). 

Millimeter and submillimeter observations of circumstellar disks can provide unique informa- 
tion. These observations probe the cool, outer parts of the disk, where giant planets are expected to 
form and contamination from the stellar photosphere is negligible. The low submillimeter opacities 
in disks can be used to extrapolate the surface density of the outer disk into the inner, optically 
thick regions and therefore determine the total disk mass (Beckwith et al. 1990). Assuming the 
submillimeter emission arises in an optically thin, isothermal portion of the disk, the flux density 
(Fi,) and disk mass (M^) are directly proportional (Hildebrand 1983): 

where d is the distance, is the opacity, and B^(Tc) is the Planck function at a characteristic 
temperature Tc. Moreover, observations and theoretical models of the opacity in the submillimeter 
indicate that Kj, is well-matched by a simple power-law in frequency with index /3, although the 
proposed normalizations vary significantly (Hildebrand 1983; Wright 1987; Pollack et al. 1994; 
Henning & Stognienko 1996). With the same optically thin, isothermal disk assumptions, the 
submillimeter continuum emission should behave roughly as oc 1/^+^. So, with major caveats 
(see §3 and the Appendix), a single submillimeter flux density can give the mass of a disk and 
> 2 flux points can reveal the frequency dependence of the opacity. Assuming a uniform grain 
composition and shape, the frequency behavior of the opacity is set by the size distribution of 
the grains in the disk. A number of single-dish surveys with single-element (or small arrays of) 
bolometers have been conducted in the Taurus-Auriga star-forming region to address these issues, 
most of which were carried out at 1.3 mm (Weintraub, Sandell, &: Duncan 1989; Beckwith et al. 1990; 
Adams, Emerson, Sz. Fuller 1990; Beckwith & Sargent 1991; Mannings & Emerson 1994; Osterloh 
& Beckwith 1995; Motte & Andre 2001). Current instrumentation provides the opportunity for 
significantly more sensitive observations of disks in the submillimeter. 

High resolution observations with (sub-)millimeter interferometers have confirmed that cir- 
cumstellar dust disks are geometrically thin with radii on the order of 100 AU (e.g., Dutrey et al. 
1996; Kitamura et al. 2002). Detailed studies of individual disks reveal molecular gas in Keplerian 
rotation around the central star (e.g., Weintraub, Masson, & Zuckcrman 1989; Koerner, Sargent, 
&; Beckwith 1993a,b; Dutrey, Guilloteau, &: Simon 1994; Koerner &: Sargent 1995; Mannings & 
Sargent 1997; Duvert et al. 1998; Guilloteau & Dutrey 1998; Simon, Dutrey, &; Guilloteau 2001; 
Corder, Eisner, &; Sargent 2005). While molecular gas is the primary reservoir of mass in a disk, it 
is difficult to directly determine from the high resolution spectral line data because the bright- 
est, easily detectable lines (i.e., the rotational transitions of CO) are optically thick (Beckwith &: 
Sargent 1993; Dutrey et al. 1996) and likely to be severely depleted (Dutrey, Guilloteau, &: Simon 
1994, 2003). Interpretation of these lines and those from trace molecular species require sophisti- 
cated models of the disk structure (e.g., Dartois, Dutrey, & Guilloteau 2003; Kamp & Dullemond 
2004) and chemistry (e.g., van Zadelhoff et al. 2001, 2003; Aikawa et al. 2002; Qi et al. 2003). 
Despite the tremendous amount of information provided by these observations, our knowledge is 
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still limited to a relatively few disks on account of the large amount of time which must be invested 
in an interferometric observation. 

Multiwavelength submillimeter data could prove useful in placing observational constraints on 
the dominant mechanism of planet formation. By comparing with infrared SEDs and diagnostics 
of accretion, we can investigate the dissipation of disks as a function of radius and see if there 
is consistency with the timescales expected from the collisional growth of planetesimals. The 
functional form of the opacity may provide information on the mean grain size distribution in the 
disk, and therefore evidence for the growth of grains demanded by the standard models of planet 
formation (e.g., Beckwith, Henning, k, Nakagawa 2000). 

In this paper, we present a large catalog of such data for most of the known YSOs in the 
Taurus- Auriga star-forming region. The survey is uniform, sensitive, and provides the most multi- 
wavelength measurements of the submillimeter continuum spectra of YSOs to date. In §2 we discuss 
the observations and data reduction procedures. In §3 we present a simple disk model and use it 
to derive circumstellar disk masses, place some new observational constraints on the submillimeter 
opacity properties of disks, and examine relationships between the disk properties and those of the 
central stars. The results are discussed in §4, and our conclusions are summarized in §5. A brief 
Appendix is included with a more in-depth discussion of the disk models we employ and comments 
on some particularly interesting sources. 

2. Observations and Data Reduction 

Simultaneous 450 and 850 /Ltm continuum photometry observations of 90 YSOs in the Taurus- 
Auriga star-forming region were obtained with the Submillimeter Common User Bolometer Array 
(SCUBA: Holland et al. 1999) at the 15 m James Clerk Maxwell Telescope (JCMT) between 2004 
February and 2005 January. Accurate reference coordinates (to ~1") for each object were obtained 
from the 2MASS Point Source Catalog. The effective FWHM beam diameters for SCUBA pho- 
tometry are 9" and 15" at 450 (Aefr = 443 /xm) and 850 /xm (Aes = 863 /xm), respectively. The 
precipitable water vapor (PWV) levels in these observations were 1.6 mm in the mean, correspond- 
ing to zenith opacities of 0.32 at 850 ;um and 1.73 at 450 /xm. More than 50% of the observations 
were conducted in very dry conditions (PWV < 1.5 mm). The data were acquired in sets of 18 s 
integrations in a small nine-point jiggle pattern with the secondary mirror chopping (typically) 60" 
in azimuth at 7.8 Hz. Each set consisted of between 15 and 40 integrations, and each source was 
usually observed for two sets. Frequent skydip observations were used to determine atmospheric 
extinction as a function of elevation and time. Pointing updates on nearby bright standard sources 
were conducted between sets of integrations: the rms pointing offsets were < 2". Mars and Uranus 
were used as primary flux calibrators, observed at least once per night when available. The sec- 
ondary calibrators HL Tau, CRL 618, and CRL 2688 were also observed approximately once every 
60 to 90 minutes. 
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The demodulated SCUBA data were fiatfielded, despiked, and corrected for extinction and 
residual sky emission using standard tasks in the SURF software package (Jenness & Lightfoot 
1998; Jenness, Lightfoot, & Holland 1998). The "unused" bolometers in the SCUBA arrays pro- 
vide a distinct advantage in sky subtraction over the standard simple demodulation utilized for 
single-element (or small array) detectors. With SCUBA, this technique has resulted in a factor 
of ~3 increase in the signal-to-noise ratio (Holland et al. 1999) and should give more robust flux 
measurements. The mean and standard deviation voltages were used to determine the flux density 
and rms noise level for each source, after appropriate scaling based on the gain values derived from 
observations of the calibrators. Repeated observations of the flux calibrators in a given night of 
observing indicate a systematic uncertainty in these gain factors of ~10% at 850 /Ltm and ~25% at 
450 fim. These systematic errors dominate the uncertainties for brighter sources. Observations of 
an additional 44 sources were obtained from the SCUBA online archive and reduced in the same 
manner, accounting for the differences in filter sets for data taken before 1999 November. In a 
few cases, observations of the same source from several different nights were combined after the 
processing to yield very sensitive data. For those cases, the combined data are consistent with the 
individual datasets when the increased integration time is considered. 

The Submillimeter High Angular Resolution Camera (SHARC-II: Dowell et al. 2003) on the 
10 m Caltech Submillimeter Observatory (CSO) telescope was also used to image 39 YSOs in the 
350 /xm continuum between 2004 March and 2005 January. The FWHM beam diameter of the 
SHARC-H point-spread function at 350 /xm is roughly 9", achieved by employing an active dish 
surface optimization system at the CSO. Due to the low atmospheric transmission at this wave- 
length, these observations were only conducted when the PWV level was < 1.6 mm, corresponding 
to 350 /um zenith opacities of less than ^^2. Opacity measurements at 225 GHz were taken every 10 
minutes with a dedicated tilting water vapor monitor observing at fixed azimuth. The observations 
were conducted by constantly sweeping the telescope in the vicinity of the source in an alt-az Lis- 
sajous pattern, providing small Nyquist-sampled maps. At least three separate maps were taken 
for each source, with between 120 and 600 s of integration per map. The aforementioned SCUBA 
calibrators were also observed every 60 to 90 minutes for pointing updates and flux calibration. The 
SHARC-H data reduction was conducted using the CRUSH software package (Kovacs et al. 2005) . 
Flux densities were measured in a circular aperture with a radius of 30", and rms noise levels were 
determined from the background pixels. Repeated measurements of standard calibration sources 
show that the absolute flux calibration is accurate to within 25%. 

Our sample was selected primarily from the compilation of Kenyon &; Hartmann (1995), and 
was designed to contain roughly equal numbers of Class H and HI objects, WTTSs and CTTSs, and 
single and multiple stars. The histograms in Figure 1 summarize some of the key properties of the 
sample. Table 1 gives a collection of submillimeter properties for 153 YSOs in Taurus-Auriga: 90 
sources with new SCUBA and SHARC-II data, 44 with archival SCUBA observations and SHARC- 
II data, 4 with data from the literature, and 15 others with SHARC-II data and flux densities from 
the literature (see the table notes). This table lists the 350, 450, 850 fim, and 1.3 mm flux densities 
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(the latter from the literature) and statistical errors (l-a rms noise levels) or 3-(T upper limits in 
units of mJy per beam, disk masses (see §3.2) and submillimctcr continuum slopes (see §3.3), and 
various other relevant properties. The projected FWHM beam diameters at the assumed distance 
of Taurus- Auriga {d = 140 pc; Elias 1978) are 1260 AU for both 350 and 450 //m and 2100 AU for 
850 fim. In this paper, we assume that all of the sources are unresolved, and therefore the values 
in Table 1 are actually the integrated continuum flux densities (in units of mJy). This assumption 
is valid for Class II and III sources, where the submillimeter emission originates in a disk with a 
radius of a few hundred AU at most (see the interferometric observations of Dutrey et al. 1996; 
Kitamura et al. 2002). 

On the other hand, submillimeter continuum maps of the Class I sources in this sample usually 
show a significant amount of extended emission from the outer envelope in addition to a bright, 
central concentration of emission (itself perhaps marginally resolved) from the disk and inner enve- 
lope (e.g.. Chandler & Richer 2000; Hogerheijde k Sandell 2000; Shirley et al. 2000; Motte & Andre 
2001; Chini et al. 2001; Young et al. 2003). The non-mapping photometry observations at 450 and 
850 /xm presented here exclude the extended emission component, and therefore only sample the 
bright peak of emission which presumably originates from warm dust in the inner envelope and/or 
a disk. Because of the unknown density structure of the inner envelope, it is not possible to unam- 
biguously determine what fraction of this emission peak is contributed by a compact object (i.e., 
disk) without interferometric observations (see the discussion by Young et al. 2003, and references 
therein). For some of these objects, there is also the possibility that the 60" chop throw would 
place the "ofF position in the extended envelope emission, and therefore the flux densities listed in 
Table 1 could be slightly underestimated. The reader should keep in mind that the submillimeter 
properties of Class I YSOs in this paper most likely refer to a combination of disk and envelope 
contributions. The notes in Table 1 provide references to submillimeter maps of the Class I YSOs 
in the literature when available. 

The primary observational goal of this survey was to take advantage of the stability and 
efficiency of the SCUBA instrument to obtain a 850 /zm sample with a relatively uniform flux 
density limit of ~10mJy (3-(t). The mean 3-a upper limit for undetected sources at 850 /um in this 
survey is 8.4 mJy (the median is the same), with a standard deviation in the upper limits of 3.1 mJy. 
For comparison, the same sources in the combined 1.3 mm surveys in Taurus- Auriga conducted by 
Beckwith et al. (1990) and Osterloh & Beckwith (1995) have a mean 3-cj upper limit of 19mJy 
(median of 16 mJy) and a standard deviation in the upper limits of 10 mJy. If we assume that 
the submillimeter continuum emission behaves as F,y oc u'^ (see §3.3), then a factor of 2.3 can be 
used to scale the 1.3 mm measurements with those at 850 /xm. The resulting scaled 1.3 mm mean 
upper limit is then 44mJy (median of 37mJy). The distributions of the upper limits of undetected 
sources are shown in Figure 2. In terms of flux density limits on undetected sources, OTir survey is 
roughly a factor of 5 more sensitive than previous single-dish work and is also considerably more 
uniform. The distributions of the signal-to-noise ratios for detected sources in the various surveys 
are similar, although there are generally higher ratios at 850 /xm. For the sources common to the 
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850 /xm and 1.3 mm samples, the detection rates are 64 ± 7% and 47 it 6%, respectively. 

3. Results 

3.1. A Simple Disk Model 

A model of the submillimeter continuum emission is needed to extract physical information 
(e.g., disk masses) from the data. In order to incorporate some non-negligible optical depth and a 
radial temperature distribution, the simplistic methods outlined in §2 (e.g.. Equation 1) are passed 
over in favor of one that fits the disk SED with a power-law structural model (see Adams, Lada, 
& Shu 1987; Beckwith et al. 1990). In this scheme, the SED (from the mid-infrared through the 
submillimeter) is assumed to be generated from thermal reprocessing of starlight by a geometrically 
thin dust disk, with the flux density given by 

cos ? f'^d 

Fu = ^J B,{Tr) (1 - e-^'^-'^^^'^') 2TTrdr (2) 

where i is the inclination angle, Tq the inner radius, Rd the outer radius, Bi,{Tr) the Planck function 
at a radius-dependent temperature, and Ty^j. the optical depth of the disk material. -"^ In essence, the 
flux density is computed by summing the thermal emission from a continuous set of dust annuli 
weighted by the radiative transfer properties of the material. The radial temperature distribution 
is taken to be a power law 

where Ti is the temperature at r = 1 AU. The optical depth is the product of the disk opacity, Ky, 
and the radial surface density profile, S^, which is also taken to be a power law: 




We assume that the opacity is a power law in frequency with index [5 and a normalization of 0.1 cm^ 
g~^ at 1000 GHz (Beckwith et al. 1990). This value assumes a 100:1 mass ratio between gas and 
dust. 

Because a given disk typically has relatively few SED datapoints, fitting the SED with the 
model described above requires that some of the remaining 8 parameters (i, Tq, Rd-, So, p, Ti, 
g, (5) be fixed. Fortunately, the precise values of i, To, and Ra do not significantly affect the 
determination of interesting physical parameters as long as they lie in a realistic range. A fiducial 



^Because Equation 2 implicitly assumes a constant source function in the disk, it is only a valid approximation 
when the inclination angle is not too large. A more sophisticated treatment of radiative transfer is required for nearly 
edge-on disks (e.g., Chiang & Goldreich 1999). 
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set of fixed parameters is adopted here: i = 0°, To = 0.01 AU, = 100 AU, and p = 1.5.^ The 
inner and outer disk radii are typical values based roughly on the dust sublimation temperature 
(e.g., DuUemond, Dominik, & Natta 2001; Muzerolle et al. 2003) and direct disk size measurements 
(e.g., Dutrey et al. 1996; Kitamura et al. 2002; Akeson et al. 2005). The surface density index, p, 
is the most difficult parameter to constrain observationally. The value selected here is obtained 
when the compositions of the planets in the solar system are augmented to cosmic abundances and 
smeared out into annuli: the Minimum Mass Solar Nebula (MMSN; Weidenschilling 1977). The 
inclination value is set merely as a computational convenience. The remaining parameters (Eo, Ti, 
q, 13) must be determined from the data. 



3.2. Disk Masses 

Submillimeter continuum observations provide measurements of disk masses. However, the 
simplistic conversion of a flux density into a mass via Equation 1 masks some important complica- 
tions. For example, could be uncertain to a factor of ~2 due to its roughly linear relationship 
with Tc- More fundamentally, the relationship between and is nonlinear due to the signifi- 
cant fraction of the submillimeter emission which is optically thick (e.g., Beckwith et al. 1990). By 
assuming optically thin emission and using Equation 1, could be underestimated (particularly 
for objects with larger flux densities). To avoid these problems and fit the SEDs with the model de- 
scribed by Equation 2, mid- and far-infrared flux densities were taken from the IRAS Point Source 
Catalog and the compilations of Weaver k, Jones (1992) and Kenyon & Hartmann (1995). Using 
data at shorter wavelengths (A < 5;um) runs the risk of contamination from an extincted pho- 
tosphere, and therefore introduces more parameters into the problem (e.g., effective temperature, 
extinction, stellar radius). The submillimeter data presented here were supplemented whenever 
possible with flux densities from the literature (Adams, Emerson, & Fuller 1990; Beckwith et al. 
1990; Beckwith & Sargent 1991; Mannings & Emerson 1994; Osterloh k Beckwith 1995; Motte & 
Andre 2001). We adopted absolute flux calibration uncertainties of 20% in the infrared and 25% 
(A < 800 jJLva) or 20% (A > 800 ixm) in the submillimeter. Systematic and statistical errors were 
combined for each individual flux density measurement. 

The disk mass and opacity index, /?, are strongly coupled parameters, making it difficult 
to independently infer their values (see Beckwith & Sargent 1991). Observations and models of 
interstellar grains in the molecular ISM, where the material is still diffuse enoTigh to safely assume 
optically thin thermal emission, indicate that /3 ~ 2 (Erickson et al. 1981; Schwartz 1982; Draine 
& Lee 1984). However, different mineralogies or grain size distributions in a disk could decrease 
the index down to /? ~ (e.g.. Pollack et al. 1994). Because of the uncertainties in independently 
measuring (3 and M^, we modeled individual SEDs for various values of (3 (between and 2), as 
well as the typical compromise value for disks, (3 = 1. Values of Ti, q, and Eo (note that for this 



^See the Appendix for a more detailed examination of the effects of various parameter choices. 
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model M^/Mq pa 2 x 10~^^So/g cm~^) for 44 objects in the sample were determined by fitting 
the SEDs to Equation 2 and minimizing the statistic. Table 2 gives the results of these fits 
for (3=1, including the reduced values (xi), degrees of freedom in the fit (v), and references 
for the infrared and submillimeter SED data from the literature. Much more sophisticated disk 
models (Men'shchikov, Henning, & Fischer 1999; Chiang et al. 2001; Semenov et al. 2005) predict 
disk masses for a few of the same sources which are within a factor of 2-3 of those presented here. 
Figure 3 shows the distributions of the best-fit values of q and Ti . The mean values of q and Ti are 
0.56±0.08 and 178±85K, respectively (quoted errors are standard deviations of the distributions). 
We define the "median disk model" to have the above set of fiducial parameters and the median 
values q = 0.58 and Ti = 148 K, as well as /3 = 1. Due to the high optical depths at infrared 
wavelengths compared to the submillimeter, the parameters of the temperature profile, Ti and q, 
are often not strongly affected by changes in /?. 

The SEDs of most of the YSOs in the survey sample were not fitted as described above 
because they either lack data (i.e., there were too few degrees of freedom), are undetected in the 
submillimeter, or have SEDs which indicate such a simple model is insufficient. However, the results 

of the SED fitting can be used to determine an empirical conversion between a submillimeter flux 
density and a disk mass. In Figure 4 we show the relationship between the 850 fim fiux densities 
and the best-fit values of (for (3 = 1) from the SED fitting. This relationship is well described 
by a simple power law, 

0.96±0.03 

(5) 



^ = (5±2)xl0-5 

M© 



(850 Aim) 



mJy 



which is shown as a solid line in Figure 4. A fit of the same data to Equation 1 (assuming the same 
opacity function given above, where Ki, = 0.035 cm^ g"-*^ at 850 /xm) is shown as a dashed line, and 
gives a best-fit characteristic temperature Tc = 20 K. For the median disk model, this value of Tc 
occurs at a disk radius of approximately 30 AU. Also shown are the relationships between Fj, and 
Md for the mean and median disk models. Disk mass values and upper limits for the objects which 
were not fitted with these models were computed from Equation 5. For sources without 850 fim 
measurements which could not be fitted with a disk model, a similar analysis as above was used to 
derive values of from the 1.3mm flux density: Md/MQ « 10~^(F,^/mJy)^-^. Disk masses (or 
3-cr upper limits) are included in Table 1 for all of the sources in the survey sample. Those masses, 
which use /? = 1, will be adopted throughout this paper, unless specifically mentioned otherwise. 
Our fitting results show that the systematic errors in the disk mass due to the a priori unknown 
value of f3 are ±0.5 dex on average, or a factor of 3, for a reasonable range of (3 (0 to 2). The 
values inferred for Class I objects should be considered only as upper limits on the disk mass, as 
there is likely a flux contribution from the inner envelope. 

Figure 5 shows the cumulative distributions of the 850 //m fiux densities and disk masses. 
This figure shows the distributions of the full sample and a subsample consisting of only those 
sources which have a > 3-cr detection at a submillimeter wavelength. The ordinates in these plots 
are defined as the probability of an object having a value equal to or greater than the abscissae. 
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in both figures, the Kaplan-Meier product hmit estimator is used to construct the cumulative 
distributions for the full sample.^ This method allows the incorporation of the 3-cr upper limits 
of the Fi, and Md values in the full sample. The computations of probabilities and their errors 
were conducted with the ASURV Rev. 1.2 software package (LaValley, Isobe, & Feigelson 1990), 
following the formalism introduced by Feigelson &; Nelson (1985). A significant caveat with these 
cumulative distributions is that there is no means to account for the uncertainties in the values of 
Fi, and M^. 

Based on the detections subsample distribution in Figure 5 (not incorporating upper limits), 
we estimate the completeness limit of the survey to be roughly 10 mJy. A log-normal distribution 
of Fi, with mean 1.20 it 0.02 (16mJy) and variance 1.08 ± 0.06 dex provides a good fit to the data 
for the full sample, whereas a mean of 1.93 it 0.01 (85mJy) and a variance of 0.41 it 0.02 dex are 
appropriate for the detections subsample. The distribution of the detections subsample in Figure 
5 indicates that 37% of the YSOs have Ma > 0.01 M©, roughly the total mass of the MMSN 
(Weidenschilling 1977). Approximately 79% of the same subsample have disks with masses greater 
than that of Jupiter. As would be expected from the relationship between Fi, and discussed 
above, the disk masses arc also log-normally distributed: the full sample with mean —3.00 it 0.02 
(10^^ Mq) and variance 1.31 it 0.06 dex, and the detections subsample with mean —2.31 ± 0.01 
(5 X 10-3 M©) and variance 0.50 ± 0.02 dex. 



3.3. Submillimeter Continuum Slopes 

The slope of the submillimeter continuum emission from a circumstellar disk is empirically 
well-described by a simple power law in frequency: Fi, oc v"'. If the emission is assumed to be 
optically thin and roughly isothermal. Equation 2 can be written Fi, oc Bi,{T)ti, oc (in the 

Rayleigh- Jeans limit). However, the submillimeter continuum has a significant contribution from 
optically thick emission originating in the dense, inner disk which causes a substantial deviation 
from the a = 2 + [3 relationship inferred for the optically thin case (Beckwith et al. 1990; Beckwith 
& Sargent 1991). We have combined the data presented here with additional flux densities from 
the literature (Adams, Emerson, & Fuller 1990; Beckwith et al. 1990; Beckwith & Sargent 1991; 
Mannings k, Emerson 1994; Osterloh & Beckwith 1995; Motte & Andre 2001) to determine the 
values of a given in Table 1. For objects with more than 2 submillimeter flux densities, a was 
measured from a linear flt in the log z^-log Fi, plane. When only 1 or 2 flux densities were available, 
values or 3-0" upper limits of a were determined from a simple spectral index. 



^Application of the Kaplan-Meier estimator to the flux density data may bo inappropriate. Because each object 
was observed either until it was detected or a rather uniform flux density limit was reached, the function which 
describes the censoring of these data is not random. Nevertheless, any effects of using the Kaplan-Meier estimator 
should only be noticed for flux densities below the completeness limit (~10mJy). Because multiwavelength SED 
data were used in determining Ma (and various SED morphologies can result in identical values of Md), the censoring 
function in that case should be randomized, and therefore the use of the Kaplan-Meier estimator is valid. 
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Figure 6 compares the submillimeter continuum slopes from 450 to 850 //m and 850 /xm to 
1.3 mm. For the full sample, the best-fit slopes are a = 2.06 it 0.02 and 1.93 ib 0.01 for the two 
wavelength regions, respectively. Because the submillimeter continuum emission is generated in the 
outer disk where temperatures are low, the Rayleigh-Jeans limit criterion is not satisfied (because 
hv ~ fcT) and the continuum slope at shorter wavelengths (nearer to the peak of the thermal 
emission) could be systematically smaller than at longer wavelengths. However, such an effect 
is not seen in Figure 6: in fact, the shorter wavelength slope is slightly steeper than at longer 
wavelengths. This implies that the shape of the submillimeter continuum is not set by the grain 
temperature distribution alone, but also by the amount of optically thick emission and/or the 
spectral behavior of the opacity function. Within the uncertainties, it does not significantly matter 
where in the submillimeter continuum the slope is measured (at least between 350 /xm and 1.3 mm). 
The slightly shallower best-fit slope for the longer wavelength data may simply be noise, or could be 
caused by a small fraction of the 1.3 mm data which are contaminated by non-disk emission: e.g., 
free- free or gyrosynclirotron radiation from a wind or outflow (e.g., Chiang, Phillips, & Lonsdale 
1996). Another possibility is a real concavity to the long- wavelength SEDs: the models of Pollack 
et al. (1994) predict a steeper opacity function shortward of ~650 ;Um. 

The cumulative distribution of a is shown in Figure 7, constructed using the Kaplan-Meier 
estimator to incorporate 3-a" upper limits. The median value of a is 2.0, while only 6% of the 
sample has a > 3, a typical value adopted in the literature due to the (incorrect) assumption of 
optically thin emission in the Rayleigh-Jeans limit with /3 = 1. A normal distribution of a with 
a mean of 1.97 ± 0.01 and a variance of 0.22 it 0.02 provides a decent fit to the data, but there is 
a slightly enhanced probability of larger a. Because some of the emission is optically thick, there 
is no straight-forward means of associating these values of the continuum slope with power law 
indices of the opacity function (/?). One approach is to allow /3 to vary in the disk SED models 
and fit it as an additional parameter (e.g., Beckwith h Sargent 1991; Mannings h Emerson 1994; 
Dutrey et al. 1996). However, in many cases this severely limits the number of degrees of freedom 
in the fits (see Table 2), which already make a number of assumptions. Beckwith & Sargent (1991) 
provide a means of relating a and /3 analytically from other parameters in the SED model fits which 
essentially indicate that (3 oc a, although there is a constant offset (see the Appendix). 

Figure 8 shows the measured values of a as a function of logM^. The shaded region on this 
diagram marks the functional form of a{M(i) for /3 = 2 which is representative of the complete range 
in the measured radial temperature distributions (see Figure 3) . The disk mass values are those for 
/3 = 1, and the error bar shown to the lower left demonstrates the systematic uncertainty introduced 
by varying (3 between and 2. All else being equal, the disk mass is roughly inversely proportional 
to the opacity, and since larger values of /? give lower opacities for a fixed frequency, a larger /? will 
also give a larger (assuming the normalization of the opacity function is fixed). Even allowing 
for such uncertainties in the disk mass and the range of temperature profiles. Figure 8 indicates 
that many of the disks in the sample have /3 < 2 (the curves showing the relationship between a 
and Md for lower values of P always fall in or below the shaded strip: see the Appendix). This 
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suggests that the large optical depths in the disk do not completely explain the shallow measured 
continuum slopes, but that there is also an evolutionary change in the typical opacity properties of 
dust grains from the ISM (where /3 = 2) to a disk. This result has also been noted in other studies 
(Beckwith & Sargent 1991; Mannings k. Emerson 1994; Dutrey et al. 1996; Dent, Matthews, & 
Ward-Thompson 1998) for various different sizes and types of samples. It is tempting to conclude 
that the apparently diminished values of /? in these disks are due to the collisional growth of dust 
grains, a necessary condition in any planet formation model (Mizuno 1980; Pollack et al. 1996). 
Models of the process indicate that grain growth can decrease [3 to values as low as zero (e.g., 
Miyake &: Nakagawa 1993; Hcnning, Michel, & Stognienko 1995). The data in Figure 8 also clearly 
show that there is no correlation between the disk mass and submillimeter continuum slope. 

The coupling of M^, /?, and the opacity normalization makes it difficult to definitively associate 
low values of the submillimeter continuum slope with decreased opacity indices (Beckwith &; Sargent 
1991). The actual value of Ki, is the main uncertainty in the conversion of a submillimeter flux 
density into a disk mass. Aside from the effect of the grain size distribution (thus the interest in 
grain growth), both the normalization and (5 depend strongly on the mineralogical composition 
of the grains (Pollack et al. 1994; Henning & Stognienko 1996) and their physical shapes (e.g., 
spherical, fractal, etc.; see Wright 1987). Further discussion of these uncertainties is given by 
Beckwith et al. (1990) and Beckwith, Henning, k, Nakagawa (2000). As an example, Wright (1987) 
indicates that can be roughly an order of magnitude higher for fractal grains compared to spheres 
at a wavelength of 1 mm (see his Figure 6) . 

Observational and theoretical uncertainties obfuscate the relationship between a measured 
submillimeter continuum slope and the functional form of the opacity in a disk. Overcoming these 
difficulties to pursue evidence of the collisional agglomeration of dust grains in the earliest stages of 
planet formation will at least require better observations, including resolved images at wavelengths 
extending beyond ~lmm (where the emission is more optically thin; e.g., Tcsti et al. 2001), fiux 
measurements near the SED turnover point (in the 100 to 300 //m range), and studies of solid-state 
dust emission features in the mid-infrared (e.g., van Boekel et al. 2004). However, our results leave 
no doubt that the measured submillimeter continuum slopes for YSOs are significantly less than 
those noted for molecular clouds, where a pa 4. 

3.4. Connections to Stellar Properties 

The physical properties of a young star and its circumstcllar disk could be related due to their 
mutual formation and subsequent gravitational and thermal links. Generally, observational indica- 
tions of any such relationship are absent, presumably due to a wide range of initial circumstcllar 
conditions for individual sources and the relatively small ranges of stellar masses and ages. As with 
Beckwith et al. (1990), we do not find any correlations between measured submillimeter properties 
and any characteristic of the stellar photosphere (e.g., effective temperature, luminosity, optical 
fluxes or colors, etc.). One perhaps notable exception is the large value of Ti derived for the two A 
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stars in this sample (AB Aur and V892 Tau). Natta, Grinin, & Mannings (2001) also suggest that 
hotter stars have generally higher dust temperatures in their disks, but there is no noticeable trend 
for the cooler (K and M) majority of this sample. This disconnect between the stellar photosphere 
and the outer disk, where the submillimeter emission is generated, is not surprising. Despite the 
increase in vertical scale height of the disk with radius expected from hydrostatic equilibrium (e.g., 
Kenyon & Hartmann 1987), radiative transfer models for structurally realistic disks indicate that 
the bulk of the submillimeter emission comes from the dust near the disk midplane, and not in 
the flared atmosphere which can be more directly affected by the stellar photosphere (Chiang k. 
Goldreich 1997, 1999). 

The gravitational link between a young star and its disk suggests that the stellar and disk 
masses may be related. Circumstellar disks are self-gravitationally stable if their mass is less than 
a fraction (a few tenths) of the stellar mass (Shu et al. 1990; Laughlin & Bodenheimer 1994). 
In principle, this could allow more massive stars to harbor more massive disks. Natta, Grinin, 
& Mannings (2001) combine interferometric measurements of disks around early-type stars with 
the 1.3 mm survey of Beckwith et al. (1990) and claim a marginal correlation between the disk 
mass and stellar mass (M*) over 2 orders of magnitude in for roughly 100 objects (however, 
see Mannings & Sargent 2000), although the dispersion is substantial. To revisit this issue, we 
have collected optical/near-infrared magnitudes and spectral classifications from various sources 
in the literature (Cohen & Kuhi 1979; Jones &; Herbig 1979; Slutskii, Stal'bovskii, & Shevchenko 
1980; Herbig, Vrba, & Rydgren 1986; Herbig & Bell 1988; Strom et al. 1989; Hartmann et al. 
1991: Gomez et al. 1992; Bouvier et al. 1993; Briccfio ct al. 1993; Hartigan, Strom, & Strom 1994; 
Martin et al. 1994; Kenyon & Hartmann 1995; Hernandez et al. 2004; White & Hillenbrand 2004, 
spectral types are listed in Table 1). A consistent set of visual extinctions was determined from 
the iy — /) color excesses, using the intrinsic colors tabulated by Kenyon & Hartmann (1995) and 
the interstellar extinction law derived by Cohen et al. (1981). De-reddened visual magnitudes and 
spectral types were converted to bolometric luminosities and effective temperatures again using 
the intrinsic values of Kenyon & Hartmann (1995). Stellar masses and ages were determined by 
reference to theoretical pre-main— sequence evolution tracks and isochrones (D'Antona &: Mazzitelli 
1997) in a Hertzsprung-Russell diagram. 

As Figure 9 demonstrates, there are no correlations between the submillimeter properties 
listed in Table 1 and stellar mass or age, but the ranges of those stellar properties (sec Figure 
1) may be too limited to infer a direct evolutionary sequence. However, the upper right panel 
of this figure shows that the region corresponding to higher mass disks at late times (> 6Myr) 
is significantly depopulated. While there are not many objects in Taurus-Auriga with such ages, 
this unoccupied region in the diagram is consistent with other studies that indicate disk fractions 
approaching zero in the 6 to 10 Myr age range (e.g., Haisch, Lada, &: Lada 2001). Figure 10 shows 
the cumulative distributions of the mass ratio of disk to star, constructed with the Kaplan-Meier 
estimator. Log-normal distributions provide poor fits in this case, but these distributions are fit 
fairly well with power laws of index between —1.5 and —2 for mass ratios larger than ~10~^. The 
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median disk to star mass ratio is 0.5%. The fraction of disks which may be self-gravitationally 
unstable (mass ratios larger than ~0.1) is negligible in Taurus- Auriga: roughly 6%, which itself 
may be an overestimate due to envelope emission for some of the Class I objects at the high end of 
the distribution. However, if /? = 2 is more appropriate, then the fraction of unstable disks can be 
as high as one third. A small fraction of objects (a few percent) has a mass ratio less than 10~^. 

3.5. The Effects of Multiplicity 

The evolution of circumstellar disks can be dictated by either internal (e.g., viscous accretion, 
gravitational instability, planet formation) or external processes. Examples of the latter include 
ultraviolet photoevaporation in the vicinity of a massive star (e.g., Johnstone, Hollenbach, & Bally 
1998), dynamical interactions with other stars in a crowded cluster environment (Kroupa 1995; 
Boffin et al. 1998) or in a local multiple star system. In the low stellar density Taurus- Auriga region, 
which is devoid of stars earlier than AO, the dominant external process affecting disk evolution is 
expected to be dynamical star-disk or disk-disk interactions in multiple star systems. Most young 
stars in nearby clusters and main-sequence stars in the field are in multiple systems, and the 
multiplicity fraction in Taurus-Auriga may be exceptionally large (Mathieu 1994; Mathicu et al. 
2000). The similar multiplicity fractions for YSOs and main-sequence field stars indicates that 
binary formation occurs early in stellar evolution (at least before the Class II stage, and likely 
much earlier), when significant circumstellar material is still present. Gravitational interactions are 
expected to severely affect the structural integrity of disks in the system, including truncation of the 
outer parts of individual circumsfe/Zar disks, gap formation in ciicumbinary disks, or even complete 
dissipation of circumstellar material via accretion or ejection (Artymowicz & Lubow 1994). 

As an example, consider a young binary system with semimajor axis a and eccentricity e, which 
also harbors two individual circumstellar disks and a larger circumbinary disk. Simulations of the 
gravitational dynamics in such a system indicate that the circumstellar disks will be truncated and 
a gap will open in the circumbinary disk, at radii which are determined primarily by the values 
of a and e (Artymowicz & Lubow 1994). The disk model described in §3.1 can be adjusted to 
determine the effects on the SED of such disk configurations by re-setting the inner and outer radii 
for the various disk components or simply by setting = for the cleared regions (e.g., Jensen, 
Mathieu, & Fuller 1996). One expected result from these SED models is that the submillimeter 
emission should be significantly diminished for systems with a projected semimajor axis (a^) on the 
order of a few tens of AU, but essentially identical to single stars for small and large Cp. Previous 
observations have indicated that < 50 — 100 AU binaries have less submillimeter emission than 
single stars or wider binaries (Jensen, Mathieu, &; Fuller 1994, 1996; Osterloh Sz Beckwith 1995), 
with the important exception of spectroscopic binaries (a < 1 AU; e.g., Mathieu et al. 1995, 1997). 

Table 3 gives a list of multiple stars in the sample and their projected separations. This 
information and the data in Table 1 have been compiled in Figure 11, which shows the 850 /xm flux 
density and disk mass as a function of projected semimajor axis. To be consistent with previous 



-15- 



work, spectroscopic binaries as well as Class I sources have been excluded in this figure and the 
analysis that follows. Unresolved higher-order multiple star systems (> 3 stars) without resolved 
observations in the literature were assigned the same Fy or Md value for all projected separations. 
Notes on assigning values for a few other systems are provided in Table 3. We utilize a variety 
of two-sample statistical tests that incorporate upper limits to determine the probabilities that 
the 850 /xm flux densities and disk masses in various binary subsamples are drawn from different 
parent distributions. The total sample is separated into categories based on projected scmimajor 
axis, resulting in three groups: close binaries with separations less than some critical value, ac\ 
wide binaries with separations larger than a^, and single stars. 

Table 4 lists the ranges of probabilities that the various subsamples for Oc = 50 AU and 100 AU 
differ from a sequence of survival analysis statistical tests performed with the ASURV software: 
the logrank, Peto & Peto, Peto &: Prentice, and Gehan tests (see the detailed descriptions by 
Feigelson & Nelson 1985). The same tests were also performed for Class II objects only. The 
results in Figure 11 and Table 4 confirm the earlier conclusions of Jensen, Mathieu, &; Fuller 
(1994, 1996): submillimeter flux densities and disk masses are significantly lower in close binaries 
(op < 50 — 100 AU) than wider or isolated systems and wide binaries essentially have the same disk 
masses as single stars. These differences arc greatest for a critical semimajor axis ac = 100 AU. 
The results for the total sample (i.e., when Class III binaries are included) generally exhibit lower 
probabilities than the subsample of only Class II objects in Table 4, with the exception of the close 
and wide binary populations with the Oc = 100 AU cutoff criterion. These differences are likely 
due to the evolutionary behavior of disks between the Class II and III stages (see §4), rather than 
an environmental effect in the multiple system. The exact probabilities for the various subsamples 
appear to be fairly sensitive to the assignment of flux densities or disk masses in unresolved higher- 
order multiple systems. High resolution interferometric observations are needed to determine the 
relative submillimeter contributions of individual components in these systems. 

In a statistical sense, the presence of a companion in the range of ~1— lOOAU decreases the 
apparent disk mass(es) in the system, presumably due to enhanced accretion onto the stars and/or 
dispersal into the local ISM. However, multiple star systems still contain disks, as evidenced by 
the relatively high detection rate in the submillimeter, 66 zt 10% for multiple systems compared 
to 58 lb 8% for single stars, as well as other inner disk signatures in the optical and infrared (e.g., 
White &: Ghez 2001). High-resolution interferometric measurements of disks in multiple systems 
have revealed a number of important exceptions to the statistical analysis above. For example, the 
GG Tau A and UZ Tau systems both have small projected separations but very large disk masses, 
the former in a circumbinary disk and the latter in a pair of disks with four stellar components 
(Koerner, Sargent, & Beckwith 1993a; Dutrey, Guilloteau, & Simon 1994; Jensen, Koerner, & 
Mathieu 1996). Moreover, single-dish continuum surveys may be missing signatures of outer disks 
in close multiple systems: a close binary (a^ = 32 AU) in the SR 24 triple system in Ophiuchus was 
found to have a large circumbinary gas disk detected in CO line emission but not in the continuum 
due to its low mass (Andrews & Williams 2005). The question of whether the large fraction of 
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young stars in multiple systems could eventually harbor planetary systems will remain unanswered 
until more detailed case studies with interferometers (e.g., Jensen & Akeson 2003) can confirm the 
properties of their disks. 

4. Discussion 

A summary of representative numbers derived from this submillimctcr survey of Taurus- Auriga 
is provided in Table 5. Listed are submillimeter detection fractions as well as median values and 
standard deviations of disk masses and submillimeter continuum slopes for the total sample and 
various subsamples of interest. Of the complete sample of 153 YSOs, 61 it 6% were detected for 
at least one submillimeter frequency, with a median ~ 5 x 10~^ and a ~ 2.0. Single 
and multiple star systems have essentially identical detection rates and similar continuum slopes. 
However, as discussed in detail in the previous section, closer binaries have statistically lower disk 
masses than wider systems or single stars. Although Figure 12 demonstrates that there is no 
direct correlation between or a and the equivalent width or luminosity of the Ha emission line, 
there is an obvious difference in the submillimeter detection fraction between WTTS and CTTS 
disks. The very high submillimeter detection rate for CTTS disks (91 it 11%) is consistent with 
all CTTSs having disk masses greater than '^10^'^ Mq. The bulk of the detected WTTS disks are 
clustered near H^(Ha) = 10 A: when the WTTS/CTTS division criterion is slightly relaxed, this 
result suggests that nearly all WTTSs are either diskless or have very low disk masses. Therefore, 
the equivalent width of the Ha emission line appears to be a fairly robust predictor of the presence 
of a "massive" disk. 

Spectral energy distribution classifications of objects in the sample were determined based on 
power-law fits from 2 to 60 /xm (when possible) with data from the literature (Strom et al. 1989; 
Weaver & Jones 1992; Kenyon & Hartmann 1995; Hartmann ct al. 2005, and references therein). 
We adopt the classification breakdown of Greene et al. (1994), using the values of the power- 
law index n (defined by uF^y oc i/") to distinguish between Class I, Flat Spectrum, Class II, and 
Class III sources. The derived classifications are listed in Table 1. Although it is not absolutely 
calibrated in time, the YSO evolution sequence defined by the shape of the infrared SED is certainly 
indicative of changes in the physical structure of the inner regions of the circumstellar disk and/or 
envelope. With the large sample of submillimeter data presented above, we can address the issue 
of corresponding changes in the physical properties of the outer disk. 

Motivated by the differences in the detection rates and median properties of the various SED 
classes listed in Table 5, the same survival analysis two-sample statistical tests used in §3.4 were 
employed to determine the probabilities that the 850 fim. flux densities, inferred disk masses, and 
submillimeter continuum slopes for various SED and Ha line strength classes are drawn from 
different parent populations. The test results are given in Table 6, and the cumulative distributions 
of Fn, Md, and a for different classes are shown in Figure 13. There are statistically significant 
progressions of decreasing submillimeter flux densities, disk masses, and continuum slopes along 
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the infrared SED evolution sequence. Flat Spectrum objects fit between Class I and II objects 

in these respects, with somewhat more similarity to the latter. Incorporating the Flat Spectrum 
objects with either the Class I or II objects does not make any significant difference in these results. 
Apparently the properties of the outer disk/envelope evolve along a similar evolutionary sequence 
as the inner disk. 

Figure 13 clearly shows that Class I objects have significantly larger submillimeter flux densi- 
ties, disk masses, and continuum slopes than Class II objects. It should again be stressed that the 
extent to which these values are representative of Class I disks, rather than disks + inner envelopes, 
is questionable. It has been suggested that many Class I disk properties could be mimicked by a 
Class II disk viewed at high inclination (e.g., Chiang &; Goldreich 1999; White &; Hillenbrand 2004). 
It is shown in Figure 14 that for a given mass, a high inclination angle produces both a lower flux 
density and a lower continuum slope; the opposite is seen in Figure 13 and Table 6. If the Class 
I emission is primarily from a disk (with only a comparatively small contribution from the inner 
envelope), a simple re-orientation of a Class II disk will not reproduce the Class I submillimeter 
properties without additional changes in mass, temperature, or opacity. The distributions found 
here of the empirical (model-independent) flux densities and continuum slopes corroborate the 
original picture of Class I sources as disk + envelope systems: the higher flux densities and "disk" 
masses may be due to additional envelope mass, and the higher continuum slopes may be due to 
the less-processed (i.e., lower amount of grain growth) dust in the envelope. A large interferometric 
sample will be required to definitively settle the issues involved in a comparison of Class I and II 
disks. 

Unfortunately, there are no measurements of a submillimeter continuum slope for any of the 
Class III objects in the sample, and the 3-a upper limits are too large to make any definitive 
statements on an evolutionary trend in a. The direct relationships between or a and the infrared 
SED slope are shown in Figure 15. There is no direct correlation with M^,^ but differences between 
the SED classes in general are apparent. A more steady decrease in a is seen across the evolution 
sequence, which if it continues would imply very shallow continuum slopes for Class III disks {a ~ 
1 — 1.5). The Spearman rank correlation coefficient in this case is —0.50, with a 99.98% confidence 
level (3.7-0"). The best-fit linear relation between the submillimeter and infrared continuum slopes 
is a = 0.40(±0.04)n + 2.09(±0.03). This trend can not be explained solely by decreasing optical 
depths in the disks along an evolutionary sequence: lower optical depths produce steeper continuum 
slopes. Another effect, such as a shallow opacity function or temperature/surface density evolution, 
must be acting to decrease a in this manner. However, interferometric observations of the Class I 
sources at several wavelengths would be required to confirm the validity of this trend. 

The submillimeter detection fraction and the fraction of objects with a near-infrared {Kg — L) 



*A significant correlation exists when L1551 IRS 5 and L1551 NE are included (the two points in the upper left 
corner), but these objects likely have large contributions to from their envelopes, and so have been excluded in 
this part of the analysis. 
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excess are identical: 60 it 7%.^ Figure 16 is a near-infrared color-color diagram that indicates the 
sources with submillimeter detections. Of the 6 sources with infrared excesses that were not detected 
in the submillimeter, 5 could have anomalous colors due to mismatched photometry and/or infrared 
companions (see the Appendix). Three of the 55 sources with essentially no near- infrared excess, 
or 5 ± 3%, were detected in the submillimeter: GM Aur, V836 Tau, and CoKu Tau/4.6 All three 
of these YSOs also have mid- and far-infrared emission, indicating that the lack of near-infrared 
excess may be due to a clearing of dust in the inner ~1AU of their disks. In addition to these 
"transition" objects, three Class III sources, V807 Tau, FW Tau, and LkHa 332/Gl, were also 
detected in the submillimeter (a 5.6 ± 3.2% detection rate), along with another possible Class III 
candidate whose SED classification remains to be confirmed due to lack of infrared data (HQ Tau) . 
However, in general a YSO with a near-infrared excess also has submillimeter emission consistent 
with a disk mass greater than ^lO^^M©, and vice versa. The small fraction of objects, less than 
10%, with evidence for an outer disk (from submillimeter data) and no inner disk suggests that the 
timescale for the disappearance of both infrared and submillimeter disk emission is relatively short; 
no more than a few hundred thousand years (i.e., < 10% of the typical YSO age in Taurus- Auriga) . 
In agreement with the comparatively low detection fraction for WTTS disks (16 =b 5%) and other 
similar analyses (e.g., Skrutskie et al. 1990; Wolk & Walter 1996; Duvert et al. 2000), these results 
imply that the inner and outer disk dissipate, or become unobservable, almost simultaneously. 

The physical mechanism responsible for the rapid and essentially radially constant "disappear- 
ance" timescale remains to be explained. Viscous accretion onto the central star alone does not 
readily produce the apparently rapid inner-outer disk dissipation (Hollcnbach, Yorke, & Johnstone 
2000). In fact, evolution under accretion processes predicts only small changes in submillimeter 
emission with time (e.g., Hartmann et al. 1998). Models which incorporate the ultraviolet photo- 
evaporation of the outer disk along with viscous accretion have more success in reproducing the 
inferred dissipation timescale, particularly for disk emission out to ~100 /xm (Clarke, Gendrin, & 
Sotomayor 2001; Armitage, Clarke, & Palla 2003). However, these "ultraviolet switch" models also 
suggest that submillimeter emission is relatively unaffected, and could therefore predict a fairly 
large fraction of WTTSs or Class III sources with submillimeter emission. Clarke, Gendrin, &: So- 
tomayor (2001) suggest that the low observed fraction of such transition objects noted by Duvert et 
al. (2000) and confirmed by the larger sample presented here may be accomodated in their models 
if different surface density profiles or viscosity values are adopted. 

An alternative explanation to actually losing disk material, onto the star or elsewhere, is a 
process which renders the dust invisible to conventional observations. A compelling possibility is 



®The values given here and in Table 5 are slightly different because a small fraction of the objects in the sample 
do not have L-band measurements in the literature. 

fourth source, BP Tau, appears to fit in this category in Figure 16. However, this is likely due to mismatched 
photometry from 2MASS and the literature, because BP Tau has a clear excess in homogeneous datasets, e.g. Kenyon 
& Hartmann (1995). 
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the collisional agglomeration of dust grains in the disk. Accelerated by gravitational settling to the 
disk midplane, the characteristic grain growth timcscales even at fairly large disk radii are thought 
to be shorter than the transition timescale inferred above (e.g., Weidenschilling & Cuzzi 1993). 
Perhaps the grain growth process has rendered the disks around many of the evolved (e.g., Class 
III) sources invisible by creating a significant population of large (~cm-sized) grains which are 
inefficient emitters at both infrared and submillimeter wavelengths. If this is to be the case, any 
collisional fragmentation process of the aggregate grains should not produce more than ~10~^ M0 
of particles which are efficient submillimeter emitters. The shallow submillimeter slopes measured 
in §3.3 and the implied low values of the opacity index /3 lend some credibility to the grain growth 
argument. Theoretical studies indicate that (3 values such as those inferred in this sample (/? ~ 1 or 
less) can be the result of a significant population of large grains (Miyake & Nakagawa 1993; Pollack 
et al. 1994; Dullemond &; Dominik 2005). The collisional growth of dust grains has also been 
inferred from submillimeter observations of both young (e.g., Beckwith & Sargent 1991; Mannings 
1994; Koerner, Chandler, & Sargent 1995) and old (e.g., Calvet et al. 2002; Hogerheijde et al. 2003) 
low mass disks, and particularly for those around the more massive Hcrbig Ac stars (Tcsti et al. 
2001, 2003; Natta et al. 2004). Complementary studies of scattered light (e.g., McCabe, Duchene, 
Sz Ghez 2003; Duchene et al. 2004b) and mid-infrared spectra (Meeus et al. 2003; Przygodda et 
al. 2003; van Boekel et al. 2004; Kessler-Silacci et al. 2005) also suggest that typical grain sizes 
are larger in these disks than for the ISM. The feasibility of this hypothesis depends critically 
on coupling with another mechanism which can diminish the accretion of gas, and therefore also 
explain the low submillimeter detection fraction for WTTSs. 

The distribution of Md shown in Figure 5 indicates that typical disks have masses significantly 
lower than those required by two of the leading theoretical models for giant planet formation. Both 
the core accretion (Pollack et al. 1996) and disk instability (Boss 1998) scenarios require disk masses 
at least a few times that of the MMSN to form a Jupiter-like planet; roughly an order of magnitude 
higher than the median mass inferred for Taurus- Auriga disks. Radial velocity surveys suggest that 
roughly 10% of stars harbor a gas giant planet within a few AU (e.g., Marcy et al. 2005), with 
the prospect that better sensitivity to long-period planets could significantly increase that fraction 
(e.g., Fischer et al. 2001). This shows that planet formation is a fairly common process. In order 
for that to be the case, the disk mass distribution constructed in §3.2 needs to be reconciled with 
the theoretical requirements of the planet formation models. Possible remedies could be extracted 
from changes to the simple disk model used in §3: for example, adjustments to the disk surface 
density profile, or a significant decrease in the normalization of the opacity function. Unfortunately, 
solutions like these will remain untested until more advanced observations become available (e.g., 
interferometers with ~0'.'l spatial resolution). A likely alternative explanation, as discussed above, 
is that a significant fraction of the disk mass is locked up in large grains or planetesimals which are 
inefficient emitters at submillimeter wavelengths. 
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5. Summary 

We have conducted a sensitive, multiwavelength submillimeter survey of 153 Taurus-Auriga 
YSOs in an effort to analyze properties of the outer regions of circumstellar dust disks. Some of 
the key results from this survey are summarized here: 

• The disk mass (or submillimeter flux density) distribution function is well matched with a 
log-normal distribution centered around 5 x 10^^ Mq with a large dispersion (0.5 dex). The 
vast majority of disks in Taurus- Auriga have substantially lower masses than is thought to be 
required for giant planet formation. However, a significant fraction of the disk mass could be 
stored in large grains or planetesimals which do not contribute to the submillimeter emission. 

• We provide the largest set of submillimeter continuum slope measurements of YSOs to date. 
The empirical behavior of the continuum from 350 jum to 1.3 mm is well-described by Fy cx 
j^2.o±o.5^ which is much flatter than for the interstellar medium. The low observed slope values 
are probably due to a combination of optical depth effects (Beckwith & Sargent 1991) and 
an inherently shallow opacity function from the top-heavy grain size distribution produced 
by collisional agglomeration of material in the disk. 

• There do not appear to be any links between stellar and disk properties in the sample, 
although the stellar masses and ages span a relatively limited range. The median disk to star 
mass ratio is ~0.5%. 

• Submillimeter flux densities and disk masses are statistically lower for stars with close com- 
panions (projected scmimajor axes less than ~100 AU) than for wider binaries or single stars. 
However, multiple star systems often still contain disks, regardless of their projected sepa- 
rations. Multiple star systems with wider separations have flux densities and disk masses 
comparable to single stars. 

• In general, the standard signatures of the inner disk (e.g., accretion diagnostics or infrared 
excess emission) arc accurate predictors of a disk mass greater than ~10~^Mq. 

• Statistically significant changes in the distribution functions of submillimeter flux densities, 
disk masses, and continuum slopes are found for the progressive stages of YSO evolution in- 
ferred from inner disk observations. These measured outer disk properties decrease from Class 
I ^ n ^ HI objects as well as for CTTSs WTTSs. The implication is that the inner and 
outer disk develop along a similar evolutionary sequence. A multiwavelength interferometric 
survey of Class I objects would be very useful for determining the relative contributions of 
a disk and inner envelope in these systems for a more sophisticated comparison with their 
presumably more evolved counterparts (Class II and III objects). 

• Only a small fraction of objects (< 10%) which have no inner disk signatures were detected in 
the submillimeter, suggesting that both infrared and submillimeter disk emission disappear on 
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a similar timescale (within ~10^ years of each other). There are two timescales in operation for 
disk evolution: (1) the relatively long (~5 to lOMyr) lifetime of Class II/CTTS disks, and (2) 
the rapid (a few x 10^ years) transition period to Class III/WTTS disks. Understanding the 
mechanisms responsible for these timescales, particularly the trigger for the transition stage, 
remain key problems in disk evolution. Some possible explanations for the essentially radially- 
independent disk dissipation timescale include viscous accretion with photoevaporation by 
the central star (e.g., Clarke, Gendrin, & Sotomayor 2001) or rapid grain growth in the early 
stages of planet formation (e.g., Weidenschilling k. Cuzzi 1993). 

We acknowledge useful conversations, suggestions, and advice from Michael Liu, Lee Hart- 
mann, Ted Simon, and Alan Boss. We would like to thank the JCMT and CSO support staffs, and 
in particular Colin Borys, for their assistance. An anonymous referee provided valuable criticism 
which helped improve this paper. This work was supported by NSF grant AST-0324328. This 
research has made use of the JCMT archive at the Canadian Astronomy Data Center, which is op- 
erated by the Dominion Astrophysical Observatory for the National Research Council of Canada's 
Herzberg Institute of Astrophysics and the NASA/IPAC Infrared Science Archive, which is oper- 
ated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 

A. Comments on Disk Models 

A simple disk structure model was presented in §3.1 and used to generate SEDs of thermal 
dust emission and derive a relationship between the submillimeter flux density and disk mass. A 
lack of BED data severely limits the number of degrees of freedom in this modeling, and therefore 
a number of parameters were fixed: inclination, inner and outer radii, and the power law indices of 
the radial surface density and frequency spectrum of the opacity. Here we determine the effects that 
reasonable alternative values for these fixed parameters would have on the Fjj — relationship. 
Readers would also benefit from the mathematical formalism discussed in detail by Beckwith et al. 
(1990). The analysis is done comparatively, relative to a fiducial disk model based on the median 
values discussed in the text: i = 0°,p= 1.5, Ti = 150 K, q = 0.6, To = 0.01 AU, Ra = 100 AU, and 
P = 1. A set of 850 /xm flux densities was computed using Equation 2 for a grid of disk masses with 
this fiducial parameter set excepting one of the previously fixed variables. The flux density — disk 
mass relationships for various values of p, q, To, Rd, Ti, and /3 are shown in Figure 17. It should 
be noted that these plots arc only intended to illustrate the effects of changing a single parameter 
in these models. In reality, the parameters are usually somewhat coupled, and therefore changes 
in one parameter affect others: such coupled effects are not considered in this simple comparative 
analysis. 

Larger values of q, To, Rd, and f3 all result in a larger disk mass for a given flux density. The 
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middle panels in Figure 17 indicate that the disk boundaries play only a minor role in setting the 
Fi, — Mfi relationship. The roles of the opacity index and temperature normalization are roughly 
those expected from the optically thin assumption given in Equation 1, where F,^ oc KyAIdT^^. 
For low-mass disks (i.e., < 10~^ M©), the shape of the radial surface density profile can have 
a significant effect on the submillimeter flux density because most of that emission is optically 
thin. For the same reason, the radial temperature behavior dictated by q has an opposite effect 
and impacts the flux densities for the more massive, optically thick disks. The largest relative 
deviations in the flux density — disk mass relationship arc set by the parameters which describe 
the radial temperature profile; Ti and q. Fortunately, the parameters of the temperature profile can 
be reliably determined for individual disks using observations in the mid- and far-infrared. Sensitive 
observations with the Spitzer Space Telescope will essentially ensure that the opacity is the single 
dominant uncertainty in the determination of a disk mass from submillimeter observations. 

In Figure 18, the same procedure as above was used to examine the effects of different fixed 
parameter choices on the relationship between the submillimeter continuum slope and the disk 
mass (note that here we have fixed /? = 1 as a reference value). Once again, the disk boundaries 
play a negligible role in the relationship. The radial surface density index has only a small impact 
on the a — relationship in general, although a constant surface density disk {p = 0) gives a 
roughly constant submillimeter continuum slope (and lower than for larger values of p) for low-mass 
disks. The parameters of the temperature profile again show the largest deviations from the fiducial 
a — relationship. The effects of a changing P in this relationship are also large. See Beckwith & 
Sargent (1991) and Mannings &; Emerson (1994) for a mathematical description of what is shown 
in these plots. 

An alternative view of the relationship between the observed continuum slope (a) and the 
power law index of the frequency behavior of the opacity (/?) can be obtained directly from these 
models. Figure 19 shows computed continuum slopes as a function of /? for various disk masses. The 
relationship is essentially linear for < 0.1 with slopes of nearly unity independent of disk 
mass, but significantly different intercept values. All of the curves fall below the nominal /3 = a — 2 
line which is representative of optically thin emission in the Rayleigh- Jeans limit. This effect is 
due to the failure of the Raylcigh-Jeans limit and the fraction of optically thick submillimeter 
emission, which can be fairly high at the shortest wavelengths (Beckwith et al. 1990; Beckwith &: 
Sargent 1991; Mannings & Emerson 1994). This synthetic grid illustrates the conclusion in §3.3 
that the roughly gaussian distribution of a centered around ~2 suggests that the opacity index is 
likely between ~0.5 and 1.5, and certainly less than the ISM value of 2 unless the disk masses are 
severely underestimated. 

Figure 20 was generated from the fiducial model (the top panel also assumes the median 
Md = 0.005 M0) described above to illustrate that the assumption of optically thin emission in the 
submillimeter continuum is not always valid (as pointed out by Beckwith et al. 1990). The disk 
becomes optically thin at the radius, ri, where t^, = Ki/E^i = 1. That criterion can be solved to 
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determine ri using Equation 4 and the relationship between So and M^i, giving 

ri = (AC,So)'/^ro = 

when p 2; inserting the parameters fixed in §3.1 gives ri ~ 184(M(^/Mq)^/"^ AU. The fraction 
of the submillimeter flux density from optically thick emission, A (defined as the ratio of the flux 
density from r < ri to the total flux density: see Beckwith et al. 1990) increases exponentially until 
A = 1 around 0.4 to 1 Jy, depending on the wavelength. Even faint submillimeter sources have 
~15% of their flux densities generated in the innermost (radially), optically thick regions of the 
disk. In terms of the fraction of the disk mass which gives rise to optically thick emission (the ratio 
of the integrated surface density from Tq to ri to the total disk mass), a more gradual trend with 
disk mass is present. Roughly 25% of the mass in a MMSN disk with the fiducial parameter set 
contributes optically thick emission in the submillimeter. Comparison of the right panel in Figure 
20 with the information in Figures 4 and 8 demonstrate the effects these relatively high optical 
depths have on estimating and a (and subsequently /3) from observations. 



{2-p)K^Md 
27r(i?2-f-r2-f) 



(Al) 



B. Comments on Individual Sources 

There are six objects in the sample with apparent near-infrared excesses, but no detected 
submillimeter emission: DP Tau, JH 223, CoKu Tau/3, FV Tau/c, CZ Tau, and CIDA-3. Three of 
these six (CoKu Tau/3, FV Tau/c, and CZ Tau) are multiple stars (see Table 3), introducing the 
possibility that the near-infrared excess is exaggerated by a red companion. The infrared colors of 
JH 223, CoKu Tau/3, CZ Tau, and CIDA-3 also seem questionable in light of the fact that these 
sources are WTTSs. Therefore, multiplicity and/or mismatched photometry from the literature for 
WTTSs could explain why submillimeter emission was not detected for most of these sources. The 
non-detection of DP Tau is more puzzling. 

As mentioned in §4, there are three sources without near-infrared excesses which are detected 
in the submillimeter (and mid- and far- infrared) : CM Aur, V836 Tau, and CoKu Tau/4. Supple- 
mentary to these transition objects arc the Class III sources detected at 850 /v,m: V807 Tau, FW 
Tau, and LkHa 332/Gl. These objects have presumably evacuated the inner portions of their disks, 
perhaps due to clearing from the planet formation process (see, e.g., Forrest et al. 2004, regarding 
CoKu Tau/4). Examinations of the SEDs for this sample highlight some other interesting sources 
in terms of possible near-infrared deficits: FM Tau, FQ Tau, UX Tau, and HK Tau. Seven late-type 
(excluding V892 Tau and SU Aur) WTTS disks are detected in the submillimeter: V773 Tau, IQ 
Tau, UX Tau, IT Tau, CoKu Tau/4, LkHa 332/Gl, and V836 Tau. A number of other sources 
detected in the submillimeter but with very little information between ~2 and 850 //m include FY 
Tau, ON Tau, CIDA-7, CIDA-8, CIDA-9, and HQ Tau. These are prime targets for mid- and 
far-infrared observations from both the ground and the Spitzer Space Telescope. One final source 
worth further investigation is CY Tau. This object has an unique SED, well-described by a shallow 
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power law from ~2/xm out to at least 1.3 mm. The models in §4 are clearly inappropriate for 
this they predict very large disk masses {Md ~ 10 Mq). Updated mid- and far-infrared 

photometry and spectroscopy of this source, along with resolved observations of the dust content 
may clear up the true nature of the circumstellar material. 
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Fig. 1. — Histograms highlighting some key properties of the submilhmeter sample. The upper left 
panel shows the distribution of spectral types from the literature (see Table 1), which is primarily 
constrained to K and M types. The upper right panel marks the relative numbers of YSOs of 
various SED classifications (Class I, Flat-Spectrum, Class II, and Class III; see also Table 1): these 
classifications are discussed in detail in §4. The bottom panels show the number distributions 
of stellar masses and ages, inferred as described in §3.4 using the D'Antona &: Mazzitelli (1997) 
theoretical pre-main— sequence models. 
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Fig. 2. — The distributions of 3-a upper limits for the sources which were not detected in both 
this survey (filled histogram) and the combined surveys of Beckwith et al. (1990) and Osterloh &: 
Beckwith (1995) at 1.3mm (unfilled histogram). The upper limits are in units of mJy, and the 
1.3 mm measurements have been scaled according to F^^ cx v"^ (see §3.3) to enable direct comparison 
with the 850 fim measurements. The submillimeter survey presented here is roughly a factor of 5 
more sensitive in terms of flux density limits, and is considerably more uniform. 
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Fig. 3. — The distributions of the best-fit disk model parameters of the normahzation (Ti: top) 
and power law index [q: bottom) of the radial temperature profile. The mean and median values 
for the sample are 0.56 it 0.08 and 0.58 for q, and 178 it 85K and 148 K for Ti (quoted errors are 
standard deviations). 
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log [Me] 

Fig. 4. — The relationship between the disk mass and 850 jim. flux density for the sources in Table 
2. The solid line is the best linear fit to the data in log-log space, given as Equation 5 in the text. 
The dashed line is the best fit to the data based on the optically thin, isothermal disk model given 
in Equation 1: the best-fit characteristic temperature is Tc = 20 K. The dotted and dash-dotted 
curves show the mean and median disk model behaviors, as described in the text. Although a 
considerable dispersion exists (the rms residual dispersion around the best-fit power law model 
given in Equation 5 is 0.2 dex), the apparent correlation here permits an empirical calibration of 
the — Fi, relationship. 



-36- 



1.0 
0.8 

^ 0.6 
_o 

^ 0.4 

0.2 
0.0 

12 3 4 

log F,{Q50 fxm) [mJy] 

1.0 

0.8 

^ 0.6 
o 

^ 0.4 

0.2 
0.0 

-5 -4 -3 -2 -1 
log [Ms] 

Fig. 5. — Cumulative distributions of 850 /xm flux densities for 140 objects in the full sample and 
the 78 of those which have > 3-a detections (top) and disk masses for 153 objects in the full 
sample and 93 objects which have submillimeter detections (bottom). The full sample distribution 
functions were computed using the Kaplan-Meier estimator to incorporate 3-cr upper limits. The 
ordinate values represent the probability that an object in the sample has a flux density or disk 
mass greater than or equal to each abscissa value. The flux densities are log-normally distributed: 
the full sample with mean 1.20 it 0.02 (16mJy) and variance 1.08 ± 0.06 dex and the detections 
subsamplc with mean 1.93±0.01 (85 mJy) and variance 0.41±0.02 dex. The same distribution holds 
for the disk masses: the full sample with mean —3.00 it 0.02 (10~^ M©) and variance 1.31 it 0.06 dex 
and the detections subsample with mean —2.31 ± 0.01 (5 x 10~^ M©) and variance 0.50 ± 0.02 dex. 
The best-fit distributions are overlaid as dashed curves. 
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Fig. 6. — The upper panels show 1.3 mm (left) and 450 //m (right) flux densities plotted against 
850 /im flux densities. The open triangles mark 3-a upper limits. For sources with no obvious error 
bars, errors are smaller than the symbol size. The solid lines are power laws of the form F,y oc v", 
and the best-fit value of a is given in the lower right corner (determined from the intercept of a 
linear fit in log-log space). The fits are only conducted for logFi,(850/xm) > 1.5, and sources with 
anomalous continuum slopes are excluded (see Table 1). The bottom panels give the residuals to the 
fits and the rms values of the residuals. Considering the scatter around the best-fit relationships, 
the data indicate that the continuum slope is roughly the same regardless of the wavelength range 
in which it is measured. 
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Fig. 7. — The cumulative distribution of submillimeter continuum slopes (a: defined by oc z/") 
for 84 YSOs, computed using the Kaplan-Meier estimator to incorporate 3-0" upper limits. The 
ordinate values represent the probability that an object in the sample has a continuum slope 
greater than or equal to a. The continuum slopes are normally distributed with mean 1.97 ± 0.01 
and variance 0.22 ± 0.02, overlaid with a dashed curve. The high slope tail of the distribution is 
slightly enhanced. 
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Fig. 8. — The submillimctcr continuum slope values {a: defined as Fi, oc v'^) plotted against the 
logarithm of the disk mass. The open triangles are 3-(T upper limits on a. A typieal error bar on 
the disk mass which incorporates the systematic uncertainties (~0.5dex) due to the unknown value 
of the opacity index is shown in the lower left corner. The dashed horizontal line marks the best-fit 
mean value of a, determined from the cumulative distribution shown in Figure 7. The dotted lines 
parallel to it mark the l-u standard deviation from the mean. The thick, solid curve marks the 
behavior of a{M(i) for the median disk model with (3 = 2. The shaded region marks the [3 = 2 
curves for extreme temperature profiles: the lower boundary for Ti = 75K and q = 0.75 and the 
upper boundary for Ti = 300 K and q = 0.35. Much of the sample appears to have (3 < 2, indicating 
a change in the opacity function in disk material relative to the interstellar medium which could 
be a signature of the collisional growth of dust grains. 



-40- 




-1.0 



1.0 



-0.5 



0.0 



0.5 



V V 



V 



• • 

•.A ^ % 



-0.5 0.0 
log M [Me] 



0.5 




-1 




1 1 

log DM97 age [Myr] 



Fig. 9. — Plots showing the relationships between disk masses or submillimeter continuum slopes 
and stellar masses or ages. Filled circles are detections and open triangles are 3-a upper limits. 
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Fig. 10. — The cumulative distributions of the disk to star mass ratio for 116 objects in the 
full sample (constructed with the Kaplan-Meier estimator) and 61 objects which were detected at 
submillimeter wavelengths. The ordinate values represent the probability that an object in the 
sample has a mass ratio greater than or equal to the abscissa value. Less than 10% of the sample 
has a mass ratio which could result in a gravitational disk instability, and even this small fraction 
may be contaminated by mass in an envelope. The median disk to star mass ratio is 0.5%. 
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Fig. 11. — The relationship between the 850 /im flux density (top) or disk mass {bottom) and 
projected semimajor axis for multiple star systems which were not resolved by the observations. 

Open triangles are S-cr upper limits, and higher order systems are connected with solid lines. The 
dashed vertical lines mark the FWHM beam radius of SCUBA at 850 /um, while the dotted lines 
mark 50 and 100 AU. For Class 11 multiple systems, the flux densities and disk masses in close 
binary systems (a^ < 50 or 100 AU) are statistically lower than those in wider binary systems or 
single stars (see Table 4). 
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Fig. 12. — The behavior of disk mass or submiUimeter continuum slope for various ranges of the 
equivalent width or luminosity of the Ha emission line. Filled circles are detections and open 
triangles are 3-0" upper limits. The dotted vertical lines mark the boundary between WTTSs and 
CTTSs. Luminosities in the Ha line were computed from equivalent widths and a linear relationship 
between the local continuum flux and the extinction-corrected i?-band flux (Reid, Hawley, & Mateo 
1995, see their Figure la). 




Fig. 13. — Evidence for outer disk evolution from the distribution functions of the 850 fim flux 
densities {top), circumstellar disk masses (middle) and submillimeter continuum slopes (bottom). 
Each cumulative distribution is labeled with the infrared SED or t^(HQ;) classification from which 

it was constructed. Statistical tests (see Table 6) confirm that there is a decrease in and 
from Class I ^ II ^ III and CTTS —>■ WTTS sources, and that a is also diminished in the evolution 
from Class I — > II sources. 
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Fig. 14. — Illustration of the effect on the relationships between the disk mass and 850 ;Um flux 
density or submillimeter continuum slope introduced by increasing the inclination angle of a fiducial 
circumstellar disk (the median disk model described in the text). Higher inclination angles (closer 
to edge-on disks) produce lower flux densities and continuum slopes for a given disk mass. 
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Fig. 15. — Relationship between the disk mass (top) or submillimeter continuum slope {bottom: 
a defined by Fi, oc z/") and the slope of the infrared SED from ~2 to 60 /xm, defined as the 
index n such that I'Fi, oc z^". The vertical dotted lines mark the borders between various SED 
classifications. Filled circles are detections and open triangles are 3-0" upper limits. There is no 
direct correlation between n and logM^. A correlation (3.7-0" with Spearman rank coefficient of 
—0.50) is seen between n and a, where the best-fit linear relationship has been overlaid as a solid 
line: a = — 0.40(=b0.04)n -|- 2.09(ib0.03). The rms residual dispersion around the best-fit line is 
0.24, which is roughly the l-o" error on a expected from absolute flux calibration uncertainties. 
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Fig. 16. — A ncar-infrarcd color-color diagram of the sample constructed from 2MASS data and L- 
band data from the literature. The solid curves mark the intrinsic colors of dwarf and giant stars. 
Dashed lines denote the reddening vectors based on the extinction law of Cohen et al. (1981), 
converted to the 2MASS photometric system. The dotted line corresponds to the classical T Tauri 
star locus derived by Meyer, Calvet, & Hilldenbrand (1997). Filled circles arc sources which are 
detected for at least 1 submillimeter wavelength (between 350 /um and 1.3 mm), and crosses are 
undetected. There are 6 objects with a near-infrared excess but no submillimeter detection which 
are addressed in the Appendix. The 4 objects with no infrared excess but a submillimeter detection 
imply that the timescales for dissipation of the inner and outer disk are similar. 
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Fig. 17. — The effects of reasonable variations of the SED model parameters on the relationship 
between the 850 /xm flux density and the disk mass. The sohd curve in each panel denotes the 
relationship for a fiducial parameter set: z = 0°, p = 1.5, q = 0.6, To = 0.01 AU, Ra = 100 AU, 
Ti = 150 K, and f3 = 1. The parameter being varied from this fiducial set is indicated in the lower 
right corner, and the different curves are marked in each panel. The shaded area is representative 
of the data values (see Figure 4). The temperature profile plays the dominant role in determining 
the relationship between Fi, and M^. 
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Fig. 18. — The effects of reasonable variations of the SED model parameters on the relationship 
between the submiUimeter continuum slope and the disk mass for a fixed opacity index P = 1. The 
solid curve in each panel denotes the relationship for a fiducial parameter set: ^ = 0°, p = 1.5, 
q = 0.6, To = 0.01 AU, Rd = 100 AU, Ti = 150 K, and (3 = 1. The parameter being varied from 
this fiducial set is indicated in the lower left corner, and the different curves are marked in each 
panel. The shaded area is representative of the data values. The lower right panel shows that /? 
is the dominant factor in setting the a — relationship, but the temperature profile also has a 
significant effect. 
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Fig. 19. — Submillimeter continuum slopes (a: defined by oc between 350 /xm and 1.3 mm 
from the fiducial disk model as a function of the power law index of the opacity for various input 
disk masses. Each curve is labeled with the logarithm of the disk mass. The dashed line marks 
the nominal relationship (5 = a — 2 for optically thin emission in the Ray leigh- Jeans limit. The 
relationship between /3 and a is linear for disk masses less than ^O.IM©, but decreased relative 
to the optically thin case due to the fraction of emission (particularly at the shortest wavelengths) 
which is optically thick. The error bar in the lower left corner shows the systematic uncertainty in 
a from absolute flux calibration errors. 
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Fig. 20. — (top): The fraction of the submillimeter flux density which is from optically thick 
emission (A) as a function of the total flux density for the flducial disk model. The various 
curves represent different wavelengths: 1.3mm (dash-dotted), 850 /xm (solid), 450 /xm (dashed), 
and 350 /xm (dotted), {bottom): The fraction of the disk mass which gives rise to optically thick 
submillimeter continuum emission at various wavelengths as a function of the total disk mass. 
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Table 1. Submillimeter Properties of Taurus- Auriga Disks' 
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Fl 
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Table 1 — Continued 
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1,3 


FM Tau 


II 


MO 


< 349 


< 442 


32 ± 8 


< 36 


0.002 


< 4.13 


3 


FO Tau 


II 


M2 




< 199 


13 ±3 


< 14 


0.0006 


< 4.29 


1,4 


FQ Tau 


II 


M2 




574 ± 170 


28 ±7 


< 40 


0.001 


4.75 ±0.61 


3,7 


FS Tau 


FS 


Ml 




< 576 


49 ±6 


< 35 


0.002 


< 3.87 


3 


FT Tau 


II 




1106 ± 82 


437 ± 56 


121 ± 5 


130 ± 14 


0.01 


1.66 ± 0.20 


3,6 


FV Tau 


II 


K5 




402 ± 90 


48 ± 5 


15 ±4 


0.001 


3.15 ±0.34 


1,4 


FV Tau/c 


II 


M4 




< 355 


< 25 


< 16 


< 0.001 




4 


FW Tau 


HI 


M4 




< 35 


4.5 ± 1.1 


< 15 


0.0002 


< 3.23 


1,3 


FX Tau 


II 


Ml 




< 169 


17 ±3 


< 30 


0.0009 


< 3.62 


3 


FY Tau 


II 


K7 




< 297 


< 27 


16 ±5 


0.0007 




3 


FZ Tau 


II 


MO 




< 273 


29 ±7 


23 ±7 


0.002 


0.55 ± 0.89 


3,7 


GG Tau 


II 


K7 


6528 ± 153 


2726 ± 250 


1255 ± 57 


593 ± 53 


0.2 


1.91 ±0.12 


1,3,6 


GH Tau 


II 


M2 




< 309 


15 ± 3 


< 30 


0.0007 


< 4.76 


1,3 


GK Tau 


II 


K7 






33 ± 7 


< 21 


0.002 




3 


GM Aur 


II 


K3 


3419 ± 133 






253 ± 12 


0.03 


2.25 ±0.23 


3,6 


GN Tau 


II 






< 187 


12 ±3 


< 50 


0.0006 


< 4.32 


1,3 


GO Tau 


II 


MO 


274 ± 26 


594 ± 185 


173 ±7 


83± 12 


0.07 


1.77 ±0.33 


3,6 


GV Tau 


I 


K3 


1676 ± 137 


1808 ± 121 


282 ± 5 


87 ±4 


0.003 


2.51 ± 0.20 


4,6 


Haro 6-13 


FS 


MO 


2729 ± 171 


1400 ± 180 


395 ± 56 


124 ± 13 


0.01 


2.15 ± 0.18 


1,3,6 


Haro 6-28 


FS 


M5 




< 2636 


11 ± 3 


< 14 


0.0006 


< 8.55 


1,4 


Haro 6-37 


II 


K6 




536 ± 204 


245 ±7 


< 88 


0.01 


1.23 ±0.62 


4 


Haro 6-39 


II 






< 903 


36 ±6 


24 ±6 


0.002 


0.95 ±0.72 


3,7 


HBC 347 


HI 


Kl 




< 138 


< 9 




< 0.0004 




4 


HBC 351 


HI 


K5 




< 166 


< 11 


< 14 


< 0.0005 




4 


HBC 352/353 


HI 


GO 




< 81 


< 9 


< 12 


< 0.0005 




4 


HBC 354/355 


HI 


K3 




< 71 


< 7 


< 16 


< 0.0004 




4 


HBC 356/357 


HI 


K2 




< 69 


< 9 


< 16 


< 0.0004 




4 


HBC 358/359 


HI 


M2 




< 72 


< 9 


< 13 


< 0.0005 




4 


HBC 360/361 


HI 


M3 




< 243 


< 14 


< 26 


< 0.0007 




1,4 


HBC 362 


HI 


M2 




< 108 


< 8 


< 14 


< 0.0004 




4 


HBC 372 


HI 


K5 




< 173 


< 8 


< 14 


< 0.0004 




4 


HBC 376 


HI 


K7 




< 49 


< 6 


< 14 


< 0.0003 




1,4 


HBC 388 


HI 


Kl 




< 54 


< 6 


< 16 


< 0.0003 




1,4 


HBC 392 


HI 


K5 




< 141 


< 6 


< 18 


< 0.0003 




1,4 


HBC 407 


HI 


G8 




< 146 


< 9 


< 64 


< 0.0004 




4 


HBC 412 


HI 


M2 




< 200 


< 9 


< 16 


< 0.0004 




4 
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Table 1 — Continued 



Object 


SED 


SpT 


A = 350 /im 


A = 450 /xm 


A = 850 


A = 1.3 mm 






a 


notes^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




(9) 


(10) 


HBC 427 


III 


K7 




< 1516 


< 14 




< 0.0007 






1 


HD 283572 


III 


G5 




< 216 


< 9 


< 15 


< 0.0004 






4 


HD 283759 


III 


F3 




< 149 


< 10 


< 35 


< 0.0005 






3 


HK Tau 


FS 


Ml 


680 ± 114 






41 ±5 


0.004 


2 


14 ±0.21 


3,6 


HL Tau 


I 


K7 


23888 ± 149 


10400 ± 1400 


2360 ± 90 


880 ± 19 


0.06 


2 


53 ±0.13 


2,3,6,8 


HN Tau 


II 


K5 




< 171 


29 ±3 


< 15 


0.0008 




< 2.79 


1,4 


HO Tau 


II 


Ml 




< 567 


44 ±6 


< 30 


0.002 




< 4.02 


3 


HP Tau 


FS 


K3 


386 ± 63 






62 ±6 


0.001 


1 


39 ± 0.24 


3 


HQ Tau 


III 






< 221 


11 ±3 


< 45 


0.0005 




< 4.72 


1,3 


Hubble 4 


III 


K7 




< 89 


< 9 


< 25 


< 0.0004 






3 


HVTau 


II 


Ml 




< 519 


47 ±4 


40 ±6 


0.002 





38 ± 0.42 


1,4,7 


IC 2087/IR 


II 






1365 ± 130 


501 ±7 




0.02 


1 


58 ± 0.42 




IP Tau 


II 


MO 




< 516 


34 ±5 


16 ±5 


0.003 


1 


80 ± 0.81 


4,7 


IQ Tau 


II 


Ml 




425 ± 26 


178 ±3 


87± 11 


0.02 


1 


53 ± 0.30 


3,6 


IS Tau 


II 


K7 




< 252 


30 ±3 


< 20 


0.001 




< 3.35 


1,3 


IT Tau 


II 


K2 




< 73 


22 ±3 


< 33 


0.002 




< 1.86 


3 


IW Tau 


III 


K7 




< 253 


< 9 


< 19 


< 0.0004 






1,4 


Jl-4423 


III 


M5 




< 52 


< 8 


< 11 


< 0.0004 






4 


Jl-4872 


III 


K7 




< 63 


< 8 


< 14 


< 0.0004 






4 


Jl-507 


III 


M4 




< 52 


< 6 


< 14 


< 0.0003 






4 


JH 56 


III 


Ml 




< 74 


< 8 


< 19 


< 0.0004 






4 


JH 108 


III 


Ml 




< 66 


< 7 


< 18 


< 0.0004 






4 


JH 112 


II 


K6 






30± 10 


< 18 


0.001 




< 5.51 


4 


JH 223 


II 


M2 




< 62 


< 7 


< 19 


< 0.0003 






4 


L1551-51 


III 


K7 




< 201 


< 13 




< 0.0006 








L1551-55 


III 


K7 




< 60 


< 5 


< 23 


< 0.0003 






1,4 


L1551 IRS5 


I 




100423 ± 812 






1276 ± 5 


0.5 


2 


95 ±0.17 


4,8 


L1551 NE 


I 




22826 ± 715 






850 ± 10 


0.3 


2 


51 ±0.24 


5 


LkCa 1 


III 


M4 




< 89 


< 8 


< 14 


< 0.0004 






4 


LkCa 3 


III 


Ml 




< 471 


< 9 


< 14 


< 0.0004 






1,4 


LkCa 4 


III 


K7 




< 37 


< 4 


< 14 


< 0.0002 






1,4 


LkCa 5 


III 


M2 




< 28 


< 4 


< 14 


< 0.0002 






1,4 


LkCa 7 


III 


K7 




< 107 


< 9 




< 0.0004 








LkCa 14 


III 


MO 




< 103 


< 9 


< 19 


< 0.0004 






4 


LkCa 15 


II 


K5 


1235 ± 80 




428 ± 11 


167 ±6 


0.05 


1 


49 ± 0.24 


4 


LkCa 19 


III 


KO 




< 90 


< 10 




0.0005 








LkCa 21 


III 


M3 




< 145 


< 10 


< 12 


< 0.0005 






4 


LkHa 332/Gl 


III 


Ml 




< 663 


12 ±3 


< 14 


0.0006 




< 6.38 


1,4 


LkHa 332/G2 


III 


K7 




< 1083 


< 9 


< 15 


< 0.0005 






1,4 


RW Aur 


II 


K3 


305 ± 32 


167 ± 34 


79 ±4 


42 ± 5 


0.004 


1 


43 ±0.21 


4 


RY Tau 


II 


Kl 


2439 ± 330 


1920 ± 160 


560 ± 30 


229 ± 17 


0.02 


1 


79 ±0.17 


2,3,6 


SAO 76411 


III 


Gl 




< 102 


< 9 


< 14 


< 0.0005 






4 


SAO 76428 


III 


F8 




< 123 


< 12 


< 14 


< 0.0006 






4 


St 34 


II 


M3 




< 243 


< 11 


< 15 


< 0.0005 






4 



Tabic 1 Continued 



Fi, [mjy] 



Object 




SpT 


A = 350 //m 


A = 450 ixm. 


A = 850 jXTSi A 


= 1.3mm 




a 


notes^ 


ft \ 


(^) 


(3) 


(4) 


(5) 


(6) 


(') 


(8) 


(9) 


(10) 


SU Aur 


II 


G2 




251 ± 40 


74 ± 3 


< 30 


0.0009 


1.73 ± 0.21 


3 


T Tau 


II 


KO 


8149 ± 253 


1655 ± 218 


628 ± 17 


280 ± 9 


0.008 


2.13 ± 0.30 


1,3,6,8 


UX Tau 


II 


K2 




523 ± 37 


173 ± 3 


63 ± 10 


0.005 


2.00 ± 0.30 


4 


UY Aur 


II 


K7 


542 ± 77 


< 523 


102 ± 6 


29 ± 6 


0.002 


2.19 ± 0.24 


4 


UZ Tau 


II 


Ml 


1823 it 142 


1811 ± 129 


560 ± 7 


172 ± 15 


0.02 


1.92 ± 0.20 


1,3,6 


V410 Tau 


III 


K3 




< 206 


7.2 ±1.8 


< 30 


< 0.0004 


< 5.27 


1,3,7 


V710 Tau 


II 


Ml 




< 495 


152 ± 6 


60 ± 7 


0.007 


2.19 ±0.36 


4,6 


V773 Tau 


11 


K3 




< 386 


9.2 ±2.9 


24 ±4 


0.0005 


< 5.88 


1,4,7 


V807 Tau 


III 


K7 




< 202 


20 ±3 


< 18 


0.001 


< 3.64 


1,4 


V819 Tau 


III 


K7 




< 317 


< 9 


< 9 


< 0.0004 




1,4 


V826 Tau 


III 


MO 




< 234 


< 7 


< 15 


< 0.0004 




1,4 


V827 Tau 


III 


K7 




< 147 


< 6 


< 19 


< 0.0003 




1,4 


V830 Tau 


III 


MO 




< 57 


< 6 


< 9 


< 0.0003 




1,4 


V836 Tau 


II 


K7 


344 ± 30 


231 ± 43 


74 ±3 


37 ±6 


0.01 


1.70 ±0.24 


4 


V892 Tau 


II 


AO 


4100 ± 560 


2570 ± 350 


638 ± 54 


234 ± 19 


0.009 


2.20 ±0.19 


2,6 


V927 Tau 


III 


M5 




< 1030 


< 10 


< 20 


< 0.0005 




1,3 


V928 Tau 


III 


Ml 




< 258 


< 8 


< 11 


< 0.0004 




1,4 


V955 Tau 


II 


MO 




< 390 


14 ±2 


< 19 


0.0005 


< 5.27 


1,4 


VYTau 


II 


MO 




< 225 


< 10 


< 17 


< 0.0005 




1,4 


Wa Tau/1 


III 


KO 




< 55 


< 6 


< 19 


< 0.0003 




1,4 


ZZ Tau 


III 


M3 




< 251 


< 8 


< 15 


< 0.0004 




1,3 



''The columns are as follows: (1) — object name; (2) — SED classification type (FS = Flat Spectrum); (3) — spectral type 
from the hterature (see §3.4); (4) — 350 ^m flux density; (5) — 450 ^m flux density; (6) — 850 fim flux density; (7) — 1.3 mm 
flux density values taken from the literature; (8) — logarithm of the disk mass (see §3.2); (9) — submiUimeter continuum 
slope (see §3.3); (10) — notes on individual sources. All flux densities are measured in units of m.Jy. Upper limits arc taken 
at the 3-(T confidence level. Quoted errors are the 1-a rms noise levels and do not include systematic errors in the absolute 
flux calibration (~25% at 350 and 450 /im, ~10% at 850 /im, and ~20% at 1.3 mm). 

''The numbers in the notes column (10) refer to the following information: 1 — 450 and 850 /im data are 
from the JCMT SCUBA archive. The original data were taken at various times between 1997 and 2002 and 
were reduced in the same way described in §2 with slight modifications for the different filter set before 1999 
November. 2 — Flux densities for V892 Tau (except 1.3 mm) and AB Aur are from Mannings (1994). Flux 
densities for HL Tau are taken from the long-term, repeated measurements at SCUBA (450 and 850 ^m: see 
http://www.jach.hawaii.edu/JCMT/continuum/calibration/sens/sccondary_2004.html and Jenness et al. 2002). Fhix den- 
sities for RY Tau (except 1.3 mm), DG Tau (except 350 ^tm), and DL Tau (except 1.3 mm) are taken from Mannings & 
Emerson (1994). The 1.3 mm flux density for BP Tau was taken from Dutrey, Guilloteau, & Simon (2003). 3 — 1.3 mm flux 
densities from Beckwith et al. (1990). 4 — 1.3mm fiux densities from Osterloh & Beckwith (1995). 5 — 1.3mm flux densities 
from Motte & Andre (2001). 6 — Additional submiUimeter flux densities (Adams, Emerson, & Fuller 1990; Beckwith & 
Sargent 1991; Mannings 1994) were used in determining a. 7 — These sources have anomalous continuum slopes, and so are 
excluded from the analysis in §3.3. Two of these objects with 850 /^m detections are known to be variable centimeter radio 
sources (V773 Tau and V410 Tau). Because V410 Tau has no signatures of disk emission throughout its SED, we assume 
that its submiUimeter flux is not from a dust disk, and instead use a 3-a upper limit of 5 mJy at 850 /xm in all the analysis. 
V773 Tau has a slight infrared excess, and so we consider the submiUimeter emission to be from the disk. The anomalous 
slopes for the other sources could be due to contamination from non-disk emission at the longest wavelengths, or errors in the 
absolute calibrations at different wavelengths. See the Appendix regarding CY Tau. 8 — References for submiUimeter maps 
of some Class I and FS YSOs in the hterature. See Young et al. (2003) for maps of 04016+2610, 04166+2706 (also Shirley 
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et al. (2000)), 04169+2702, 04239+2436 (also Chini ct al. (2001)), 04248+2612, 04295+2251, 04302+2247, 04361+2547, and 
04381+2540 (also Hogerheijde & Sandell (2000)). Sec Hogcrhcijde & SandcU (2000) for maps of 04325+2402 and 04368+2557 
(also Chini ct al. (2001)). Chandler & Richer (2000) provide maps of HL Tau and Sandell & Weintraub (2001) provide maps 
of L1551 IRS 5. Maps of the extended submillimeter emission around T Tau are provided by Weintraub et al. (1999). 
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Table 2. Results of SED Fits^ 



Object 




1 


Ma [Mq] 




A 


r,2 

Xv 


u 




notes 


(1) 






(4) 






(5) 


(o) 


lis 




(9) 


AA Tau 


129 


0.56 


1.3 ±0.2 X 


10- 


-2 


0.31 


4.7 


8 


15 


1,2,3 


AB Aur 


367 


0.45 


4.4 ±0.6 X 


10" 


-3 


0.17 


2.6 


6 


65 


1 2 


BP Tau 


117 


0.64 


1.8 ±0.4 X 


10" 


-2 


0.42 


2.0 


4 


11 


1 2 


CI Tau 


152 


0.56 


2.8 ±0.6 X 


10- 


-2 


0.36 


1.7 


9 


15 


12 3 


CW Tau 


204 


0.62 


2.4 ±0.4 X 


10- 


-3 


0.24 


5.7 


6 


27 


2,3 


DE Tau 


130 


0.55 


5.2 ± 0.8 X 


10- 


-3 


0.25 


0.6 


3 


19 


1,2 


DF Tau 


126 


0.74 


4± 1 X 


10- 


-4 


0.28 


3.6 


3 


23 


1,2 


DG Tau 


288 


0.51 


2.4 ±0.3 X 


10" 


-2 


0.29 


2.7 


13 


40 


12 3 4 


DH Tau 


109 


0.55 


3.3 ±0.7 X 


10- 


-3 


0.24 


2.4 


4 


19 


1 2 


DK Tau 


175 


0.70 


5± 1 X 


10- 


-3 


0.34 


0.6 


4 


18 


1 2 


DL Tau 


149 


0.62 


9±2 X 


10- 


-2 


0.58 


0.6 


12 


10 


2,3,4 


DM Tau 


111 


0.51 


2.4 ±0.4 X 


10- 


-2 


0.34 


2.0 


6 


14 


2,3 


DN Tau 


117 


0.60 


2.9 ± 0.6 X 


10- 


-2 


0.43 


1.5 


5 


11 


2,3 


DO Tau 


193 


0.52 


7± 1 X 


10- 


-3 


0.23 


2.3 


12 


33 


1,2,3,4 


DQ Tau 


143 


0.60 


1.9 ±0.5 X 


10- 


-2 


0.36 


0.4 


4 


14 


2 3 


DR Tau 


216 


0.58 


1.9 ±0.3 X 


10- 


-2 


0.32 


1.5 


10 


27 


12 3 4 


DS Tau 


97 


0.67 


6±1 X 


10- 


-3 


0.39 


1.0 


4 


11 


1 2 


FT Tau 


121 


0.58 


1.4 ±0.2 X 


10- 


-2 


0.34 


5.0 


8 


13 


2,3 


FV Tau 


190 


0.53 


1.1 ± 0.2 X 


10- 


-3 


0.17 


1.4 


4 


38 


2 


FX Tau 


118 


0.63 


9±2 X 


10- 


-4 


0.24 


0.6 


2 


21 


1,2 


FZ Tau 


148 


0.73 


2.0 ± 0.6 X 


10- 


-3 


0.33 


1.9 


2 


17 


2 


GO Tau 


172 


0.56 


2.3 ±0.8 X 


10- 


-1 


0.76 


1.0 


8 


10 


2 3 


GM Aur 


136 


0.44 


2.5 ±0.5 X 


10- 


-2 


0.28 


2.1 


6 


30 


2 3 


GO Tau 


90 


0.62 


7±2 X 


10- 


-2 


0.62 


1.2 


7 


7 


2 3 


GV Tau 


339 


0.46 


2.8 ±0.3 X 


10- 


-3 


0.16 


1.5 


7 


78 


2,3 


Haro 6-13 


181 


0.47 


1.1 ± 0.1 X 


10- 


-2 


0.23 


0.9 


11 


33 


2,3,4 


HK Tau 


140 


0.44 


4.5 ±0.5 X 


10- 


-3 


0.18 


2.8 


9 


28 


12 3 5 


HL Tau 


277 


0.42 


6.5 ±0.8 X 


10- 


-2 


0.36 


0.7 


12 


33 


12 3 5 


HN Tau 


164 


0.56 


8±2 X 


10- 


-4 


0.18 


1.9 


2 


33 


1,2 


HP Tau 


201 


0.47 


1.0 ± 0.3 X 


IQ- 


-3 


0.14 


5.9 


4 


160 


1 2 


IP Tau 


107 


0.63 


2.8 ± 0.6 X 


10- 


-3 


0.29 


1.1 


3 


15 


2 


IQ Tau 


121 


0.60 


2.2 ± 0.3 X 


10- 


-2 


0.40 


2.4 


5 


12 


1,2,3 


IT Tau 


104 


0.61 


1.5 ± 0.3 X 


10- 


-3 


0.25 


5.1 


2 


17 


1,2 


LkCa 15 


117 


0.52 


4.8 ±0.9 X 


10- 


-2 


0.41 


2.2 


5 


13 


1,2 


RYTau 


342 


0.66 


1.8 ±0.3 X 


10- 


-2 


0.34 


1.4 


13 


52 


1,2,3,4 


SU Aur 


264 


0.48 


9±3 X 


10- 


-4 


0.14 


3.2 


4 


66 


1,2 


T Tau 


338 


0.45 


8.2 ±0.9 X 


10- 


-3 


0.19 


1.9 


12 


62 


1,2,3,5 


UX Tau 


132 


0.41 


5.1 ± 0.7 X 


10- 


3 


0.17 


3.3 


5 


31 


1,2 


UY Aur 


226 


0.53 


1.8 ± 0.3 X 


10- 


-3 


0.18 


1.8 


5 


47 


1,2 


UZ Tau 


167 


0.61 


6± 1 X 


10- 


-2 


0.49 


2.1 


9 


13 


1,2,3 


V710 Tau 


112 


0.58 


1.9 ±0.4 X 


10- 


-2 


0.37 


0.6 


5 


12 


2,3 


V836 Tau 


97 


0.58 


1.0 ±0.3 X 


10- 


-2 


0.34 


0.6 


5 


12 


2 


V892 Tau 


461 


0.58 


9±1 X 


10- 


-3 


0.25 


1.3 


8 


58 


2,3 


V955 Tau 


139 


0.59 


5± 1 X 


10- 


-4 


0.19 


6.2 


2 


27 


2 



^The table columns are: (1) — object name; (2) — best-fit value of the temperature 
at 1 AU in K (typical errors are ± a few K); (3) — best-fit value of the radial power law 

index of the temperature profile (typical errors are ± 0.02); (4) — best-fit value of the 
disk mass for /3 = 1 ; (5) — fraction of the 850 yaa flux density from optically thick regions 
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in the disk; (6) — reduced chi-squarcd statistic; (7) — number of degrees of freedom in 
the fit (i.e., number of datapoints — number of fitted parameters [= 3]); (8) — inferred 
characteristic temperature in K from inversion of Equation 1; (9) — notes on literature 
sources for the data as follows; 1 = 10 jum photometry from Kenyon & Hartmann (1995), 
2 = IRAS photometry (12, 25, 60, 100 /im) from Weaver & Jones (1992), 3 = submil- 
limeter photometry from Beckwith & Sargent (1991), 4 = submillimeter photometry from 
Mannings & Emerson (1994), 5 = submillimeter photometry from Adams, Emerson, & 
Fuller (1990). The largest values of xi usually are due to structure in the infrared SED, 
where the errors on flux densities are low, or slightly inconsistent absolute calibration in 
the submillimeter. 
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Table 3. Multiple Star Systems^ 



Object 


Sp ["] 


rcfb 


ap [AU] 


Object 


6p ["] 




ap [AU] 


(1) 


(2) 


(3) 


{^) 


(1) 


(2) 


(3) 


(4) 


r\A 1 1 Q 1 O'TV^Q 


4.U 


1 Q 
lb 


ODU 


11 rv lau 


O A 

Z.4 


( 


o4U 


H/iQoc; 1 Ofino 


8.2 


18 


1140 


HN Tau 


3.1 


8 


430 


r\A1AQ -i-ltn o 

U4tz4i5-r z o±Z 




1 Q 
1<5 


ft/in 
D4U 


rHJ iau 


ft o 


1 1 
11 


y ( u 


r^^T^-.i T^ni-. /Q 

L^oiYU lau/o 


O 1 

z.l 


( 


zyU 


rlr^ iau 


U.Uz 


1 (\ 

lU 


Q 


L^Zj iau 


U.oo 


Q 
O 


A ft 

4D 


H V lau 


n HQ A n 
U.Uo, 4.U 


6 


4, 560 


DD Tau 


0.56 


8 


/ y 


io lau 


0.22 


9 


31 


DF Tau 


n no 
U.Uy 


"1 ft 

Id 


lo 


i 1 1 au 


O A 

z.4 


•J 
I 


o4U 


iJix iau 


Z.O 


D 


OOU 


IW lau 


0.28 


8 


39 


L)y^ iau 


SD 


1 A 
14 


U.Uo 


J l-4c5 / Z 


Q Q 
O.O 


1 ft 

ID 


4DU 


Uo lau 


V 1 

( . 1 


1 1 
11 


yyu 


T 1 c;i I'D G r; 
LilOOl llto O 


U.oU 


1 


/10 

4z 


FF Tau 


n HQ 
U.Uo 


lU 


/I 
4 


ijK(^a o 


SD, U.4o 


11, 16 


0.03, 67 


r (J lau 


U. lo 


Q 
O 




ijKL^a / 


1 n 
l.U 


Q 


1 A{\ 

14U 


r lau 


n 70 
U. / o 


Q 

o 


IIU 


LiKllQ; ooz/ IjI 


0.22 


8 


31 


r o iau 


U.zo 


Q 
O 


OO 


LiKilQ: OOz/ijrZ 


U.ztt 


Q 




r V lau 


U. < z 


Q 
O 


lUl 


HW Aur 


1 A 

1.4 


I 


onn 
zUU 


r V lau/c 


U. /4 


C 

D 


1U4 


ot o4 


SD 


on 
zU 




r w lau 


n 1 o Q 
U.lD, Z.o 


1^, Id 


OO Qon 


1 iau 


n 1 n n "rn 
U.IU, U. ( U 


17 Q 
1 ( , 


1 A nc 
14, yo 


FX Tau 


0.90 


9 


130 


UX Tau'^ 


0.14, 2.7, 5.9 


15, 7, 7 


20, 380, 830 


GG Tau'* 


0.26 


8 


36 


UY Aur 


0.88 


9 


120 


GH Tau 


0.33 


9 


46 


UZ Tau'^ 


sb, 0.35, 3.7 


19, 9, 7 


0.1, 49, 520 


GK Tau 


2.5 


7 


340 


V410 Tau 


0.07, 0.29 


16 


10, 41 


GN Tau'* 


0.04 


6 


6 


V710 Tau 


3.2 


7 


450 


GV Tau 


1.2 


3 


170 


V773 Tau 


0.11 


9 


15 


Haro 6-28 


0.66 


8 


92 


V807 Tau 


0.02, 0.37 


12, 9 


3, 52 


Haro 6-37 


2.7 


8 


370 


V826 Tau 


sb 





0.05 


HBC 351 


0.61 


8 


85 


V892 Tau 


4.0 


13 


560 


HBC 352/3 


8.6 


8 


1200 


V927 Tau 


0.29 


8 


41 


HBC 354/5 


6.3 


11 


880 


V928 Tau 


0.18 


8 


25 


HBC 356/7 


2.0 


11 


280 


V955 Tau 


0.33 


8 


46 


HBC 360/1 


7.2 


8 


1010 


VY Tau 


0.66 


8 


92 


HBC 412 


0.70 


8 


98 


ZZ Tau 


0.03 


12 


4 



*The table columns are: (1) - object name; (2) - projected separation in arcseconds for each pair (sb 
denotes a spectroscopic binary); (3) - reference for each separation measurement (sec fe); (4) - projected 
semimajor axis in AU, assuming a distance of 140 pc (spectroscopic binary separations are the best-fit 
values of asini). 

''Projected separations, dp, are average values from the following sources: = Mundt ct al. (1983); 1 = 
Rodriguez et al. (1986); 2 = Simon et al. (1987); 3 = Leiiiert & Haas (1989); 4 = Haas, Leinert, & Zinnecker 
(1990); 5 = Leinert et al. (1991); 6 = Simon et al. (1992); 7 = Reipurth & Zinnecker (1993); 8 = Leinert et 
al. (1993); 9 = Chez, Neugebauer, & Matthews (1993); 10 = Richichi et al. (1994); 11 = Mathieu (1994); 
12 = Simon et al. (1995); 13 = Leinert, Richichi, & Haas (1997); 14 = Mathieu et al. (1997); 15 = Duchene 
et al. (1999); 16 = White & Chez (2001); 17 = Tamazian (2004); 18 = Duchene et al. (2004a); 19 = Martin 
et al. (2005); 20 = White & Hillenbrand (2005). 

■^These high-order multiple systems were assigned various flux densities and disk masses for the analysis 
described in §3.5 and represented in Figure 11 and Table 4. T Tau — The stellar components are too close to 
rule out a large circum-triple disk from current interferometric observations, so each separation is assigned 
the same flux/mass. UX Tau — High-resolution 1.3 mm observations show that A is the only component 
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with a disk (Jensen & Akeson 2003), so the 2 separations (ap = 380 and 830 AU) are assigned the same 
flux/mass value, and the close binary B is given a flux density equal to the completeness limit (lOmJy) 
and a corresponding disk mass from Equation 5. UZ Tau — We adopt the result of Jensen, Koerner, & 
Mathieu (1996) that the W close binary contributes ~19% of the total flux density in the system, and scale 
the flux densities and disk masses accordingly for each separation. HV Tau — We assume that C is the 
sole component with a disk because the AB close binary has a Class III SED, and C has been shown in the 
optical to have an edge-on disk (Stapelfeldt et al. 2003). The AB separation {ap = 4 AU) is assigned a flux 
density equal to the completeness limit (10 mjy) and a corresponding disk mass from Equation 5. 

"^GG Tau; a quadruple source, these numbers refer only to the A close binary, as the B system falls outside 
the 850 /xm SCUBA beam. GN Tau: White & Chez (2001) note a significantly larger projected separation, 

Op = c.'aa. 
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Table 4. Multiplicity Effects on Disk Properties^ 



samples 




Oc [AU] 






close vs. wide 




50 


85 


-91% 


33 - 54% 


close vs. single 




50 


9 


-60% 


< 7% 


wide vs. single 




50 


53 


-76% 


53 - 80% 


close vs. wide 




100 


97 


-99% 


> 99% 


close vs. single 




100 


26 


-70% 


34 - 83% 


wide vs. single 




100 


77 


-91% 


65 - 89% 


Class II, close vs. 


wide 


50 


90 


-94% 


42 - 73% 


Class II, close vs. 


single 


50 


75 


-90% 


16 - 33% 


Class II, wide vs. 


single 


50 


4 


- 14% 


89 - 96% 


Class II, close vs. 


wide 


100 


94 


-98% 


> 98% 


Class II, close vs. 


single 


100 


83 


-90% 


91 - 98% 


Class II, wide vs. 


single 


100 


12 


-38% 


18 - 38% 



''Results of censored statistical tests to determine how stel- 
lar companions affect circumstellar disks. See also a graphical 
comparison in Figure 11. Two-sample tests were performed on 
the categories listed in the left column and described in the 
text (§3.5): close binaries (with < Oc), wide binaries (with 
Op > flc), and single stars. The probabilities that the 850 jum 
flux densities (P(F„)) or disk masses (P(M<j)) are drawn from 
different parent populations are given in the last two columns. 
The ranges in the probabilities are representative of the various 
statistical tests. 
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Table 5. Summary of Submillimeter Properties^ 



sample 




/smm 




median [Mq] 


o-(Md) [dex] 


No, 


median a 


a{a) 


Class I 


16 


1.00 ±0.25 


16 


3 X 10-2 


0.59 


16 


2.51 


0.54 


Flat-Spectrum 


9 


1.00 ±0.33 


9 


4 X 10-3 


0.53 


7 


2.02 


0.70 


Class II 


74 


0.86 ± 0.11 


64 


3 X 10-3 


0.67 


34 


1.79 


0.51 


Class III 


54 


0.07 ±0.04 


4 


6 X 10-4 


0.26 









WTTS 


61 


0.15 ±0.05 


9 


2 X 10-3 


0.65 


5 


1.73 


0.27 


CTTS 


74 


0.91 ±0.11 


67 


4 X 10-3 


0.72 


40 


2.02 


0.62 


multiples 


61 


0.66 ±0.10 


40 


2 X 10-3 


0.74 


19 


2.13 


0.52 


singles 


92 


0.58 ±0.08 


53 


1 X 10-2 


0.64 


38 


2.01 


0.60 


total 


153 


0.61 ±0.06 


93 


5 X 10-3 


0.71 


57 


2.01 


0.57 



'^The N values are the total numbers of sources in each subsample. /smm is the fraction of those objects 
which were detected at a submillimeter wavelength along with the 1-<t Poisson counting error. Median values 
of disk masses and submillimeter continuum slopes axe given, along with the standard deviations (it) of those 
values. The a{M^) values are on a log scale. 



Table 6. Outer Disk Evolution: Survival Analysis Tests^ 



samples 






P(a) 


Class I vs. Flat-Spectrum 


63 - 85% 


> 99.2% 


74 - 87% 


Class I vs. Class 11 


> 99.96% 


> 99.9999% 


> 99.7% 


Class I vs. Class III 


> 99.9999% 


> 99.9999% 




Flat-Spectrum vs. Class II 


89 - 94% 


60 - 83% 


1 - 36% 


Flat-Spcctrum vs. Class III 


> 99.9999% 


> 99.9999% 




Class II vs. Class III 


> 99.9999% 


> 99.9999% 




WTTS vs. CTTS 


> 99.9999% 


> 99.9999% 


26 - 65% 



''The results of survival analysis tests as described in the text (§4). The 
values of P{Fy), P(M^), and P(ci) arc the probabilities that the 850 /im 
flux densities, disk masses, and submillimeter continuum slopes, respec- 
tively, of the two samples in the first column are drawn from different 
parent populations. The ranges of P values are representative of the var- 
ious statistical tests. See Figure 13 for a graphical representation of this 
comparison. 



